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PREFACE 


This  report  describes  the  work  performed  under  Naval  Surface  Warfare  Center 
(NSWC)  Contract  No.  N00167-95-C-0033,  “NDE  of  Ceramic  Matrix  Composites  for  Ra- 
domes”,  SBIR  Phase  II.  Analytical  studies  reported  in  Part  1  were  performed  by  Materials 
Sciences  Corporation  (MSC)  staff;  Dr.  Sailendra  N.  Chatterjee  (principal  investigator).  Dr. 
Chian-Fong  Yen  and  Mr.  Dominic  Mirto.  Experimental  studies  reported  in  Part  2  were  per¬ 
formed  under  the  supervision  of  Professor  John  C.  Duke,  Jr.,  Department  of  Engineering 
Science  and  Mechanics,  Virginia  Polytechnic  Institute  and  State  University  (VPI&SU).  Dr. 
John  M.  Liu  was  the  technical  coordinator  at  NSWC. 
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INTRODUCTION 


The  goal  of  this  effort  was  to  choose  one  or  more  nondestructive  techniques,  which  are 
simple  and  yet  can  provide  quick  evaluation  of  next  generation  radomes  for  quality  verifica¬ 
tion  during  process  development  or  for  product  acceptance  purposes.  Reasonable  quantita¬ 
tive  predictions  of  physical  property  variations  over  the  entire  part  should  also  be  possible 
from  test  data  or  generated  images.  It  may  be  noted  that  currently  used  inspection  tech¬ 
niques  using  transmitted  light  are  not  suitable  for  reinforced  ceramics  which  do  not  allow 
such  transmission.  Phase  I  studies,  summarized  in  the  next  section,  indicate  that  optical, 
infrared,  dielectric  or  ultrasonic  (thru-transmission)  imaging  can  be  useful  for  determining 
changes  in  void  or  reinforcement  content  commonly  encountered  in  reinforced  ceramics. 
Large  variations  of  void  and/or  inclusion  content  appear  to  be  the  dominant  type  of  anomaly 
in  unreinforced  as  well  as  reinforced  ceramics  and  can  significantly  influence  the  electro¬ 
magnetic  performance  of  radomes.  X-ray  radiography  may  be  useful  for  identifying  other 
kinds  of  defects. 

Studies  were  conducted  in  this  phase  to  test  the  usefulness  of  various  inspection  tech¬ 
niques  for  particle  reinforced  ceramics  being  considered  for  next  generation  radomes. 
Works  performed  include  development  of  measurement  techniques,  procedures,  scanning 
devices  and  softwares  for  nondestructive  imaging  of  dome  shaped  radomes  using  ultrasound 
and/or  microwaves.  Such  inspections  can  be  employed  at  critical  steps  of  manufacturing 
for  quality  assurance  and  particularly  prior  to  system  assembly  and  performance  testing  as 
specified  by  the  Navy  CDRL  and  other  acceptance  documents.  Various  analyses  were  also 
performed  to  evaluate  the  performance  of  electromagnetic  sensors  used  for  nondestructive 
imaging  and  to  estimate  the  effects  of  patch  defects  on  radome  performance. 
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BACKGROUND 


Radome  materials  proposed  for  future  generations  of  tactical  missiles  will  require  stabil¬ 
ity  at  high  temperatures  (to  1400°C),  low  dielectric  constants  and  low  loss  tangents.  Also, 
they  will  need  dimensional  stability  (low-thermal  expansion  coefficients),  moderately  good 
strength  and  toughness  (to  withstand  handling,  aerodynamic  and  thermal  stresses)  as  well 
as  thermal  shock  and  erosion  resistance  (Ref.  1,  2).  Moisture  absorption  should  be  minimal. 
Finally  the  density  should  be  as  low  as  possible  to  reduce  weight  and  increase  payloads. 
Continuous  fiber  (in  laminated  or  woven  forms)  as  well  as  particle  or  whisker  reinforced  ce¬ 
ramics  (Ref.  3,  4)  are  being  considered  for  such  applications,  since  these  composites  can  be 
tailored  to  provide  desired  dielectric  and  thermo-mechanical  properties.  Since  fabrication 
processes  for  ceramic  matrix  composites  are  new  and  complicated,  it  is  possible  that  distri¬ 
bution  of  reinforcements  and  other  microstructural  details  in  a  fabricated  part  or  product 
may  not  be  the  same  as  desired.  In  addition,  such  components  may  also  contain  flaws  like 
surface  cracks,  voids,  breaks  in  continuous  fibers,  reinforcement/matrix  separations,  etc. 
Nondestructive  evaluations  are  useful  for  locating  and  quantifying  these  features  as  well  as 
their  effects  on  performance  and  the  expected  life  of  a  product.  They  can  also  serve  a  role 
in  the  development  of  improved  processes  and  fabrication  techniques. 

Phase  I  studies  were  aimed  at  classifying  possible  types  of  flaws  and  anomalies  in 
whisker  reinforced  ceramics  being  considered  for  radome  applications,  obtaining  some  ana¬ 
lytical  or  semi-empirical  estimates  of  the  effects  as  well  as  critical  sizes  and  locations  of 
such  defects  and  developing  some  nondestructive  testing  techniques  which  may  be  em¬ 
ployed  for  their  characterization.  Although  it  is  advantageous  to  determine  as  many  proper¬ 
ties  (or  anomalies)  as  possible  using  the  same  technique,  it  is  known  from  NDE  studies  on 
ceramics  (Ref.  5-7)  that  sometimes  multiple  techniques  may  be  needed  to  quantify  certain 
types  of  anomalies.  Various  NDT  methods  were  employed  for  examination  of  whisker  rein¬ 
forced  ceramic  samples.  A  brief  description  of  the  results  of  Phase  I  studies  (Ref.  8)  fol¬ 
lows. 

Void  content  and  changes  in  amount  of  whiskers  were  found  to  be  the  most  important 
anomalies  in  whisker  reinforced  ceramics.  Crack  like  flaws  can  also  occur  in  such  materials. 
Following  remarks  can  be  made  about  their  influence  on  material  properties  and  characteri¬ 
zation  by  NDE. 
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1  Analytical  results  indicate  that  26%  porosity  can  cause  about  42%  reduction  in  Young's 
and  shear  moduli  of  whisker  reinforced  ceramics  (Figure  1).  Such  reductions  are  possibly 
acceptable.  Also,  11%  porosity  can  cause  greater  than  17%  reduction  in  modulus, 
which  should  be  detectable  by  ultrasonics  with  reasonable  accuracy.  In  addition,  analyti¬ 
cal  results  also  indicate  that  a  10%  increase  in  fiber  content  causes  8%  increase  in 
moduli,  which  should  be  detectable.  However,  change  in  moduli  due  to  void  or  fiber  con¬ 
tent  may  offset  each  other,  which  may  pose  a  problem  for  NDT. 


C* =100  GPa,  G«r=38  GPa,  E^=220  'GPa,  <w=S7  GPa 


Figure  1.  Elastic  Moduli  of  Mullite/Celsian  With  Voids  Th.:  Plane 
Isotropic  Through  the  Thickness  Properties 


2.  Thermal  conductivity  and  loss  tangent  are  affected  by  fiber  and  void  content  and  those 
changes  may  not  be  of  concern.  Changes  in  conductivity  due  to  porosity  can  be  quite 
significant  (Figure  2)  and  hence  infrared  scanning  may  be  an  effective  NDT  tool. 

3.  Dielectric  constant  increases  slightly  with  fiber  content,  but  it  is  more  strongly  affected 
by  porosity  (Figure  3).  About  12%  reduction  occurs  due  to  10%  increase  in  void  con¬ 
tent.  It  is  not  known  whether  such  a  change  is  acceptable  from  the  point  of  view  of  ra- 
dome  performance  (effects  on  signals  due  to  variation  in  void  content).  However,  NDE 
techniques  like  eddy  current  scanning  with  a  ferromagnetic  backing  material  which  em¬ 
ploy  changes  in  dielectric  properties  of  the  ceramic  composite  as  the  key  parameter,  may 
be  able  to  detect  this  change  with  reasonable  accuracy. 
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Figure  2.  Thermal  Conductivity  of  Mullite/Celsian  With  Voids: 
Plane  Isotropic  Properties 
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Figure  3.  Dielectric  Constant  of  Mullite/Celsian  With  Voids 
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4.  Strength  is  strongly  influenced  by  void  content  (Figure  4),  but  may  not  be  sensitive  to 
fiber  content.  Isolated  crack  like  defects  can  cause  additional  reduction  in  strength 
(Table  1).  However,  sizes  of  such  defects  to  cause  strength  reduction  are  comparatively 
large.  For  example,  in  a  material  with  1 1  %  void  content,  crack  sizes  have  to  be  greater 
than  0.15  mm  to  reduce  the  strength  below  the  level  of  about  80  MPa.  Nondestructive 
techniques,  such  as  radiography  should  be  able  to  detect  such  flaws  with  reasonable  ac¬ 
curacy. 


Figure  4.  Effect  of  Porosity  on  Bending  Strength  of  Celsian  Containing 
Classified  and  Unclassified  Mullite  Whiskers 

Table  1 .  Flaw  Characteristics  to  be  Detected  in  Mullite/Celsian 


Allowable  Strength, 

MPa 

Maximum  Void  Content 

Percent 

a.  Crack  of  Depth  * 

Greater  than,  mm 

80 

11.3 

0.15  -  0.3  mm 

60 

17.7 

0.25  -  0.5  mm 

40 

26.7 

0.6  -1.2  mm 

*  In  materials  of  thickness  greater  than  5  time  this  depth 
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5.  Infrared  scanning  and  optical  transmission  appear  to  have  good  potential  for  develop¬ 
ment  of  practical  methods  for  efficient  nondestructive  evaluation  of  samples  as  well  as 
dome  shaped  products.  However,  quantitative  data  for  full  or  subscale  components  are 
required  before  they  can  be  considered  for  practical  application. 

6.  Eddy  current  scanning  with  ferro-magnetic  backing  plates  may  also  be  useful. 

7.  X-ray  radiography  should  be  useful  for  quantifying  crack  like  defects,  although  no  tests 
were  conducted  in  Phase  1  for  this  purpose.  Use  of  radiography  as  attempted  in  Phase  I 
appears  to  offer  limited  potential  due  to  difficulties  either  with  practical  implementation 
or  data  interpretation. 

8.  Ultrasonic  thru-transmission  scans  indicate  changes  in  wave  velocities  caused  by  fiber 
and  void  content  differences  (Figure  5).  Use  of  ultrasonic  or  other  types  'of  scanning 
should  be  useful  for  dome  shaped  parts  (Figure  6). 


Figure  5-  Thru-Transmission  Ultrasonic  Scan  of  40%  Mullite  Whiskers 
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From  the  results  summarized  above,  it  is  clear  that  significant  changes  occur  in  sound 
wave  velocities,  thermal  conductivities  as  well  as  dielectric  constants  due  to  variations  in 
void  content.  Also,  nondestructive  examinations  conducted  in  Phase  I  show  that  ultrasonic, 
infrared  and  dielectric  as  well  as  optical  imaging  can  provide  quick  quantitative  evaluation  of 
reinforced  ceramic  samples  for  determining  variations  in  related  material  properties  due  to 
voids.  Crack  like  flaws  or  large  voids  can  have  significant  influence  on  strength  and  other 
NDT  techniques,  such  as  radiography,  are  needed  for  characterizing  such  flaws.  Further  de¬ 
velopment  of  one  or  more  promising  methods,  especially  for  quantitative  mapping  of  prop¬ 
erty  variations  over  an  entire  component  was  the  goal  of  the  Phase  II  work  described  in  the 
following  sections  and  in  Part  2. 
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PHASE  II  TECHNICAL  OBJECTIVES 


Results  of  Phase  I  studies  (Ref.  8)  summarized  in  the  previous  section  indicate  that  opti¬ 
cal,  infrared,  dielectric  or  ultrasonic  (thru-transmission)  imaging  can  be  useful  for  determin¬ 
ing  changes  in  void  or  reinforcement  content  commonly  encountered  in  reinforced  ceramics. 
Presence  of  voids  appears  to  be  the  dominant  type  of  anomaly  in  unreinforced  as  well  as 
reinforced  ceramics  (Figure  7).  X-ray  radiography  may  also  be  useful  for  identifying  different 
kinds  of  anomalies.  However,  usefulness  of  such  techniques  have  to  be  demonstrated  be¬ 
fore  they  can  be  used  for  routine  application  for  inspection  of  reinforced  ceramics  being 
considered  by  the  Navy.  The  specific  objectives  of  the  Phase  II  effort  as  described  in  the 
original  proposal  are  outlined  below. 


Figure  7.  Voids  and  Porosities  in  Subscale  Celsian  radomes  (X-Ray  Tomogram) 
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(i)  To  procure  reinforced  ceramic  samples  and  components  which  are  being  considered  for 
Navy  applications. 

<ii)  To  obtain  analytical  estimates  of  effects  of  particle  reinforcements,  voids  and  cracks  on 
various  material  properties  such  as  stiffness,  strength  and  dielectric  properties  and  criti¬ 
cal  flaw  characteristics  as  well  as  the  influence  of  dielectric  property  changes  on  ra- 
dome  performance. 

(iii)  To  conduct  analytical  studies  of  some  NDE  methods  such  as  dielectric  imaging  for  use 
and  interpretation  of  NDE  data. 

(iv)  To  perform  NDE  of  samples  using  optical,  infrared,  dielectric  and  ultrasonic  techniques 
and/or  other  methods  as  needed. 

(v)  To  conduct  NDE  of  components  using  the  most  promising  methods  for  detailed  imaging 
over  the  entire  component  with  appropriately  selected  data  acquisition  techniques. 

(vi)  To  perform  some  destructive  tests  to  generate  material  property  data. 

(vii)  To  correlate  data  from  nondestructive  and  destructive  tests  with  analyses. 

(viii)To  deliver  representative  hardware  and  software  for  the  most  promising  NDE  methods 
and  to  suggest  a  plan  of  inspection. 

(ix)  To  suggest  procedures  for  the  development  of  accept/reject  criteria. 
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PHASE  H  TASKS 


The  following  is  a  list  of  tasks  originally  proposed  for  the  Phase  II  effort. 

1 .  Selection  of  materials  and  procurement  of  samples/components. 

2.  Obtain  estimates  of  reinforcements,  voids,  cracks  or  other  flaws/anomalies, 
their  effects  on  changes  in  material  properties  as  well  as  radome  performance. 

3.  Analytical  studies  in  support  of  NDE. 

4.  NDE  of  samples. 

5.  Data  acquisition  during  NDE  of  components. 

6.  Limited  destructive  tests. 

7.  Correlation  of  destructive  tests  and  NDE  data. 

8.  Inspection  plan  and  representative  hardware. 

9.  Suggestions  for  accept/reject  criteria. 

Tasks  6  and  7  were  originally  proposed  for  a  whisker  reinforced  ceramic  (Mullite-Celsian) 
for  which  test  material  was  readily  available.  By  the  time  the  program  was  started  new  ma¬ 
terials  under  development  (Zr-PSBN)  were  selected  for  investigation.  Some  samples  were 
manufactured  for  NDE  studies,  but  the  number  was  not  enough  to  generate  a  reasonable 
data  base.  This  part  of  the  report  (Part  1 )  summarizes  the  works  performed  in  the  program 
and  the  details  of  work  and  relevant  results  from  Tasks  1 ,  2,  3,  and  9.  Details  of  sample 
preparation,  experimental  studies,  inspection  plan  and  hardware  (Tasks  1,  4,  5,  and  8,  per¬ 
formed  under  the  supervision  of  Professor  John  C.  Drake  of  VPI&SU)  are  described  in  Part  2 
of  the  report. 
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DESCRIPTION  OF  WORK 


TASK  1.  MATERIALS  AND  REINFORCEMENTS 

One  material  which  is  being  considered  for  next  generation  radomes  is  a  Phosphate 
Bonded  Silicon  Nitride  developed  by  the  ceramics  group  at  NSWC,  White  Oak.  Ease  in 
processing,  good  water  stability,  low  and  thermally  stable  dielectric  constant  and  loss  tan¬ 
gent  are  some  of  its  attractive  properties.  Raw  materials  required  are  Zirconia  powder 
(30%)  and  Silicon  Nitride  particles  (70%).  Phosphoric  acid  is  used  as  a  bonding  agent 
which  reacts  with  Zirconia  to  yield  some  Zirconium  Phosphate.  It  may  be  noted  that  the 
final  material  may  be  viewed  as  Silicon  Nitride  Particles  bonded  by  Zirconium  Phosphate, 
Residual  Zirconia  and  voids.  Residual  Zirconia  has  a  very  high  dielectric  constant,  but  the 
resulting  material  has  a  large  percentage  of  voids  which  makes  the  effective  dielectric  con¬ 
stant  of  the  material  quite  low.  Some  samples  (with  and  without  flaws)  and  small  compo¬ 
nents  (domes,  shown  in  Figure  8)  were  manufactured  by  Jonathan  Medding  (VPI&SU)  at 
NSWC  facilities  and  were  interrogated  at  VPI  laboratories  (see  Part  2  for  details).  Scale  up 
is  being  attempted  by  Lockheed-Martin  Corporation  (Dallas,  TX)  under  a  contract  awarded 
by  NSWC. 

Other  materials,  such  as  Pyroceram  and  slip  cast  fused  silica  are  also  being  considered 
for  next  generation  radomes. 
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Figure  8.  Small  Domes  Manufactured  Compared  with  a  Full  Size 
Dome  of  a  Different  Material 
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TASK  2.  ESTIMATES  OF  PROPERTY  VARIATIONS  AND  EFFECTS  ON  RADOME 
PERFORMANCE 

Subtask  2.1  Effects  of  Voids  and  Inclusions  on  Material  Properties 

The  base  material  (Phosphate  Bonded  Silicon  Nitride)  investigated  also  contains  about 
20%  voids  and  5%  residual  Zirconia.  Presence  of  voids  is  desirable,  since  it  reduces  the  di¬ 
electric  constant.  On  the  other  hand,  residual  Zirconia  has  a  very  high  dielectric  constant 
and  any  increase  in  its  volume  percentage  will  increase  the  dielectric  constant  of  the  final 
material.  In  addition,  changes  in  void  and/or  Zirconia  content  over  the  entire  radome  may 
influence  its  performance.  For  these  reasons,  some  estimates  of  the  influence  'of  void  and 
Zirconia  contents  are  of  interest.  Such  estimates  can  be  obtained  using  closed  form  analyti¬ 
cal  results  for  a  two  phase  medium  (Ref.  9),  where  the  medium  is  treated  as  a  composite 
sphere  assemblage.  For  the  purpose  of  evaluation  of  the  changes  caused  by  voids,  it  is  as¬ 
sumed  that  the  voids  are  spherical  and  are  distributed  over  the  volume  of  a  parent  material. 
Since  the  dielectric  constant  of  Phosphate  Bonded  Silicon  Nitride  with  residual  Zirconia  is 
not  known,  it  was  estimated  assuming  that  a  typical  material  has  a  dielectric  constant  of 
5.0  (measured  on  actual  samples  by  the  researchers  of  the  Ceramic  Group  at  NSWC,  Ref. 
10).  Also,  because  there  may  be  some  uncertainty  in  the  actual  void  content,  the  dielectric 
constant  of  the  parent  material  was  estimated  for  three  values  of  void  volume  (20,  25  and 
30  percent)  and  these  estimates  were  employed  to  find  expected  changes  in  the  dielectric 
constant  with  void  volume.  The  results  are  shown  in  Figure  9. 

Using  the  same  approach,  effects  of  Zirconia  content  on  the  dielectric  constant  were  es¬ 
timated  for  three  values  of  the  base  material  (with  voids)  property  and  the  results  are  shown 
in  Figure  10. 

Based  on  the  results  shown  in  Figure  10,  the  dielectric  constant  of  Phosphate  Bonded 
Silicon  Nitride  can  be  estimated  and  the  result  is  given  in  Table  2. 
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Effect  Of  Spherical  Voids  On  Phosphate  Bonded  SiaN4 


Estimates  Based  On  Composite  Sphere  Assemblage,  5%  Change 
!n  Voids  Causes  7%  Change  In  Dielectric  Constant 


Figure  9.  Changes  in  Dielectric  Constant  with  Void  Volume  for  Three  Assumed  Values 
of  Void  Content  in  a  Typical  Phosphate  Bonded  Silicon  Nitride 


Effect  Of  Spherical  Residual  ZrO,  Particles 

Estlmotes  Based  On  Composite  Sphere  Assemblage,  5%  Change 
In  ZrOi  Volume  Causes  11%  Change  In  Olelecfrlc  Constant 


Figure  10.  Changes  in  Dielectric  Constant  with  Residual  Zirconia  Content  for  Three 
Values  of  Base  Material  (With  Voids)  Property 
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Table  2.  Dielectric  Constants  of  Constituents 


Material 

Dielectric  Constants 

Silicon  Nitride/Zirconium  Phosphate 

5.8 

Voids  (Air) 

1.0 

Zirconia 

40.0 

Typical  zirconium  phosphate  bonded  si 

icon  nitride  contains  about  20%  voids 

and  5%  residual  Zirconia.  Such  a  material  has  a  dielectric  constant  of  about  5.0. 

Variations  of  dielectric  constant  as  well  as  other  properties,  such  as  modulus  of  elastic¬ 
ity,  are  also  useful  for  judging  whether  NDE  techniques  can  be  used  effectively  for  deter¬ 
mining  variations  in  void  and/or  Zirconia  content.  Modulus  of  elasticity  of  a  two  phase  solid 
can  be  evaluated  using  the  methods  employed  in  Phase  I  studies.  Table  3  summarizes  the 
results  of  changes  in  material  property  for  a  5%  change  in  void  or  Zirconia  content.  The 
strength  variation  reported  in  Table  3  is  based  on  an  empirical  relation  obtained  in  Phase  I 
for  a  Mullite  reinforced  Celsian  material.  The  same  relation  was  used  to  find  the  effect  of 
varying  void  content  on  strength  reported  in  Table  4.  Table  4  also  lists  the  sizes  of  surface 
cracks,  which  must  be  exceeded,  to  cause  additional  strength  reduction  based  on  an  as¬ 
sumed  value  of  fracture  toughness  (kIC  =  1 .8  MPa  4rn  ).  Based  on  these  results,  it  is  pos¬ 
sible  to  identify  variation  in  void  content  and  minimum  surface  crack  depths  which  should 
be  detected  by  NDE.  These  values  are  underlined  in  Table  4. 

The  results  of  NDE  of  discs  with  defects  performed  in  the  program  (see  the  summary  of 
works  performed  in  Task  4  in  a  later  section  and  in  Part  2  of  the  report)  show  that  variations 
of  porosity  can  indeed  be  quantified  by  ultrasonic  imaging.  Using  the  results  discussed  ear¬ 
lier,  the  changes  in  the  dielectric  constant  in  the  porous  region  can  be  estimated.  Micro- 
wave  NDE,  using  a  flanged  co-axial  line  or  a  rectangular  wave  guide  can  also  yield  an  esti¬ 
mate  of  the  changes  in  the  dielectric  constant.  Such  changes  can  have  a  significant  influ¬ 
ence  on  the  electromagnetic  performance  of  the  radome,  such  as  boresight  errors  or  receiv¬ 
ing  patterns  of  the  antenna.  The  next  subtask  describes  the  works  performed  in  an  attempt 
to  estimate  these  influences. 
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Table  3.  Summary  of  Results  for  a  Typical  Phosphate  Bonded  Silicon  Nitride 


5%  Increase  in  Void  Volume  Fraction 

7%  Decrease  in  Dielectric  Constant 

10%  Decrease  in  Young’s  Modulus 

5%  Decrease  in  Density 

20%  Decrease  in  Tensile  Strength 

5%  Increase  in  Residual  Zirconia 

Volume  Fraction 

1 1  %  Increase  in  Dielectric  Constant 

May  also  influence  modulus,  density,  and  sonic 

velocity 

•  Some  changes  in  thermal  conductivity  and  thermal  expansion  are  also 
expected 


Table  4.  Some  Variations  and  Flaws  to  be  Detected 


Void  Content 

Change  in 

Void  Volume 

Strength 

(MPa) 

Minimum  Surface 
Depth  in  mm  to 
Cause  Additional 
Strength  Reduction* 

20%  + 

0 

88 

0.16-0.37 

22.1% 

2.1 

80 

0.20-0.45 

25.1% 

5.1 

70 

0.25  -  0.58 

28.5% 

8.5 

60 

0.35  -  0.80 

*  Based  on  assumed  fracture  toughness  1 .8  MPa  -yfm  for  all  cases. 

+  Assumed  as  representative  material 

Underlined  variations/flaws  should  be  detectable  bv  NDE 


Subtask  2.2  Effects  of  Patch  Defects  on  Radome  Performance 


The  primary  function  of  the  radome  is  to  provide  an  enclosure  for  the  antenna.  It  would 
be  ideal  to  have  a  radome  with  sufficient  strength  to  withstand  environmental  effects  and 
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aerodynamic  loads,  but  with  a  low  dielectric  constant  to  allow  electromagnetic  waves  to 
pass  through  it  with  very  little  changes  in  direction  and  attenuation.  However,  since  the 
dielectric  constant  of  the  radome  materials  is  always  greater  than  that  of  air  (unity),  the  ap¬ 
parent  direction  of  arrival  of  a  target  return  can  be  different  from  the  true  direction.  The 
bias  of  the  apparent  direction  towards  azimuth  (xa)  (or  elevation  (yA))  axes  of  the  antenna  is 
called  the  azimuth  (or  elevation)  boresight  error,  which  is  usually  measured  in  milliradians 
(MRAD).  Boresight  error  varies  with  the  orientation  of  the  direction  of  target  return  (zA  axis 
of  the  antenna)  with  respect  to  the  radome.  Estimates  of  boresight  error  variation  and  its 
slope  (with  respect  to  the  change  in  the  direction  or  gimbal  angle)  based  on  calculations 
and/or  experimental  measurements  are  often  used  to  determine  the  correct  direction  of  the 
target  return.  As  discussed  earlier  variations  of  dielectric  properties  are  expected  in  ceramic 
or  ceramic  composite  radomes  over  a  finite  sized  region  (a  patch)  due  to  changes  in  porosity 
or  reinforcement  content.  Such  variations  will  obviously  influence  the  boresight  error  and  it 
is  of  interest  to  obtain  an  estimate  of  these  variations  due  to  such  a  patch  type  defect 
(Figure  11).  Experimental  measurements  with  patches  of  different  sizes  and  different  dielec¬ 
tric  properties  would  involve  an  elaborate  program  and  were,  therefore,  beyond  the  scope  of 
this  effort.  Therefore,  in  this  subtask  an  attempt  was  made  to  estimate  the  influences  of 
patch  type  defects  using  the  computer  code  RACP™  (Radome  Analysis  Computer  Program, 
Ref.  11)  obtained  from  scientific  Software,  Maitland,  Florida. 

The  RACP™  code  is  quite  versatile  and  allows  evaluation  of  boresight  errors  for  different 
types  of  radomes  with  variable  thickness  and  properties  as  well  as  monopulse  antenna  near 
fields  (receiving  patterns)  for  various  polarizations  of  the  incident  plane  wave  and  the  an¬ 
tenna.  However,  in  its  present  form  the  code  does  not  allow  for  consideration  of  a  patch 
type  defect  in  the  radome.  The  code  was,  therefore,  modified  to  consider  such  property 
variations  over  a  region  bounded  by  two  prescribed  axial  coordinates  (z)  and  circumferential 
coordinates  (<t>)  as  shown  in  Figure  12.  Since  it  is  expected  that  the  property  change  will  be 
maximum  at  the  center  of  the  patch  and  will  reduce  to  zero  at  its  boundaries,  it  is  assumed 
that  the  property  change  (real  and  imaginary  part  of  the  dielectric  constant  or  thickness)  is 
of  the  following  form 


A  =  Ao  cos 


-cos 


0 


(1) 
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Figure  1 1 .  Radome  with  a  Patch  Defect,  p.  8x 


Figure  1 2.  Dielectric  Property  Validation  and  Position  of  the  Patch  p.  8y 


where  zi  and  <j>i  are  the  coordinates  of  the  center  of  the  patch  and  l  and  0  quantify  its  di¬ 
mension  in  axial  and  circumferential  direction.  It  may  be  noted  that  <j>i  defines  the  relative 
angular  position  of  the  patch  with  respect  to  the  antenna.  4>i  =  90°  when  an  axial  plane 
(through  radome  z-axis)  passing  through  the  center  of  the  patch  coincides  with  the  azimuth 
plane  (xAzA)  of  the  antenna  (Figure  12).  For  =  0°,  the  axial  plane  coincides  with  yAzA  or 
elevation  plane  of  the  antenna  (Figure  1 2). 

A  typical  tangent  ogive  radome  shown  in  Figure  1 1  was  considered  for  demonstration 
purposes.  The  following  geometric  and  material  properties  were  used  (Figure  13). 

Radome  base  diameter,  Dos  =  1 2.88  inch 

Radome  length,  Los  =  26.43  inch 

Bulkhead  z  coordinate,  zBot  =  0.0  inch 

Antenna  Aperture  Radius,  R*p  =  5.0  inch 

Gimbal  point  from  radome  base,  R®  =  7.057  inch 

Antenna  aperture  offset  from  gimbal,  RA  =  1 .908  inch 

Frequency,  F  =  1 1 .803  GHZ 

Radome  wall  thickness,  DIN  =  0.351  inch 

Dielectric  constant  er,ER  =5.0 

Loss  tangent,  TD  =  0.03 

For  estimation  of  changes  in  boresight  errors  due  to  the  patch,  the  dielectric  constant,  sIr  at 
center  was  assumed  to  be  4.0  (Figure  12).  For  estimation  of  changes  in  receiving  patterns 
one  more  value  of  sr  at  the  center  of  the  patch  (sr  =  2.0)  was  also  selected. 

Patches  of  trapezoidal,  but  nearly  square,  shapes  (dimensions  in  the  z  direction  equal  to 
that  in  circumferential  direction  at  <j)  =  <j>i)  were  considered.  Figure  12  illustrates  such  a 
patch  of  dimension  6"  in  both  directions.  In  the  results  presented  here  and  in  Appendix  A, 
where  additional  results  are  given,  a  defect  of  size  6"  (or  4")  implies  such  a  nearly  square 
shaped  patch. 
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Figure  13.  Geometry  of  Dimensions  of  Radome  and  Antenna 
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Figure  1 4  shows  the  variation  of  the  azimuth  boresight  error  for  a  vertically  polarized  an¬ 
tenna  and  incident  wave  with  yaw  gimbal  angle  0l  (the  xAzA  plane  of  the  antenna  lies  in  the 
plane  of  the  paper  in  Figure  13)  and  the  patch  size.  It  may  be  seen  that,  as  expected,  for 
positive  values  of  0l  (normal  to  antenna  plane  zA  moving  towards  the  defect)  the  azimuth 
boresight  error  increases  with  the  patch  size.  A  patch  size  higher  than  6"  causes  the  bore- 
sight  error  to  increase  tenfold  (for  0L  *  8°),  but  below  6"  it  increases  about  three  times  and 
may  be  acceptable  if  appropriate  corrections  are  made  to  account  for  the  presence  of  the 
defect.  On  the  other  hand.  Figure  14  shows  that,  for  obvious  reasons,  the  boresight  errors 
are  less  influenced  for  negative  values  of  0L. 


Figure  14.  Effect  of  Defect  Size  at  90  Degree  Position  and  z  =  15  in.  on  Bore 
Sight  Azimuth  Error  for  Vertically  Polarized  Antenna  and  Incident  Wave 


Figure  1 5  shows  the  variation  of  the  azimuth  boresight  error  for  a  right  hand  circularly 
polarized  antenna  and  a  vertically  polarized  incident  wave  for  a  fixed  patch  of  size  6"  (and 
for  no  defect)  with  yaw  gimbal  angle  0L  and  different  angular  position  of  the  defect 
(<j>i  =  90°  to  0°).  Although  the  boresight  errors  change  quite  a  bit  for  all  values  <f>i,  but  they 
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may  still  be  acceptable.  However,  the  final  decision  regarding  acceptability  can  only  be 
made  after  some  tests  are  performed  with  some  patches  of  different  sites  and  locations  and 
the  data  are  correlated  with  the  analytical  results. 


— *-NoD*hd 
■  ■  ■  900  d*cr»e,  6  Inch  D«ftct 
— Ar-67.5  tftgrw.1 6 Inch  D«f*ct 
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— 0.0  tfcgrx.  t  Inch  Otfact 


Figure  15.  Effect  of  Angular  Position  of  6  in.  Defect  at  z  =  15  in.  on  Bore  Sight 
Azimuth  Error  for  RHC  Polarized  Antenna  and  Vertically  Polarized  Incident  Wave 


Effects  of  a  6"  size  patch  with  er  =  4.0  and  2.0  on  the  receiving  pattern  of  the  antenna 
(sum  channel  amplitude  and  phase  at  azimuth  cut)  are  shown  in  Figures  16  and  17.  Some 
changes  are  noticeable,  but  the  differences  from  no  defect  case  for  er  =  4.0  are  less  than 
that  for  s,  =  2.0. 

Additional  results  and  the  details  of  changes  in  the  input  data  required  to  consider  the 
effects  of  patch  defects  are  described  in  Appendix  A. 
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Figure  16.  Effect  of  Defect  on  Receiving  Pattern  -  Sum  Channel  Amplitude  Azimuth  Cut 


— 4 — 6“  Defect,  «p«  too  -  2 
—A —  €'  Defect,  ef'olon  "  4 
— •—No  Defect _ _ 


Figure  1 7.  Effect  of  Defect  on  Receiving  Pattern  -  Sum  Channel  Phase  Azimuth  Cut 
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TASK  3.  ANALYTICAL  STUDIES  OF  NDE  METHODS 

Since  the  primary  interest  is  in  the  electromagnetic  properties  of  the  material,  it  should 
be  useful  to  determine  variation  of  the  dielectric  constant  (or  other  related  responses)  of  the 
material  over  an  entire  part  using  NDE  techniques.  Open-ended  wave  guide  type  sensors  are 
attractive  devices  suggested  by  various  investigators.  In  this  technique,  the  material  (a  ce¬ 
ramic  plate  or  shell)  is  interrogated  using  a  flanged  open-ended  rectangular  wave  guide  or  a 
co-axial  line  and  the  changes  in  the  actual  or  normalized  admittance  of  the  wave  guide  (or 
voltage  standing  wave  ratio,  VSWR  and  phase  of  the  reflection  coefficient)  are  measured. 
These  changes  can  be  related  to  those  in  the  dielectric  constant  of  the  material  due  to  voids 
(and/or  delaminations).  There  are  many  other  factors  which  affect  the  admittances,  such  as 
backing  material  (on  the  other  side  of  the  wave  guide)  and  any  gap  between  the  flange  and 
the  material.  Analytical  studies  and  results  are  useful  to  determine  the  effects  of  such  pa¬ 
rameters,  such  that  suitable  care  can  be  taken  during  measurement.  In  addition,  the  results 
can  be  used  for  analysis  and  interpretation  of  data.  Therefore,  a  computer  code  has  been 
written  to  perform  such  analytical  calculations.  Mathematical  procedures  and  numerical 
methods  employed  are  presented  in  Appendix  B  of  the  report. 

Appendix  B  also  describes  the  input  data  required  for  approximate  evaluation  of  the  ad¬ 
mittance  of  the  rectangular  wave  guide  or  the  co-axial  probe  used  for  interrogation  of  a  lay¬ 
ered  dielectric  medium.  A  listing  of  the  computer  code  is  also  included.  Figures  18  and  19 
show  the  changes  in  the  normalized  conductance  and  susceptance  of  rectangular  wave 
guide  sensor  (of  dimensions  0.9"  x  0.4")  due  to  the  expected  variations  in  the  dielectric 
constant  of  a  typical  Phosphate  Bonded  Silcon  Nitride  plate  of  thickness  0.2"  backed  by  air. 
Results  are  given  for  two  cases  (no  air  gap  and  an  air  gap  of  0.02"  between  the  sensor  and 
the  plate)  and  three  frequencies.  The  results  show  that  the  charges  in  susceptance  are  very 
sensitive  to  the  dielectric  constant  at  a  frequency  of  8  GHz  (Figure  18)  and  are  also  strongly 
influenced  by  the  presence  of  the  air  gap.  Therefore,  a  frequency  in  the  neighborhood  of  8 
GHz  may  be  suitable  for  interrogating  the  material  and  careful  consideration  should  be  given 
to  the  presence  of  unavoidable  air  gap  during  measurements.  The  effect  of  the  air  gap  on 
measurements  using  a  co-axial  probe  and  correlation  of  such  test  data  with  analytical  results 
are  discussed  later. 
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Normalized  Conductance 


(  Air  Backed  Plate  Of  Thickness  0.2",  fan  6  =  -.003  ) 


a  **  8.0  GHz.  Airgap  •  0“ 

<u  m  fl.O  GHz,  Ajrgop  •  0.02“ 

«  *■  10.0  CHz,  Atrgop  *  0" 
o  —  10,0  GHz.  Airgap  •  0.02" 
«  «  12.0  GHz,  AJrgap  *  0" 
a  m  12.0  GHz,  Airgap  •«  0.02" 


Dielectric  Constant,  c 


Figure  18.  Normalized  Susceptance  vs.  Dielectric  Constant  of 
Phosphate  Bonded  Silicon  Nitride 


(  Air  Backed  Plate  Of  Thickness  0.2",  fan  6  =  -.003  ) 


Dielectric  Constant,  t 

Figure  19.  Normalized  Conductance  vs.  Dielectric  Constant  of 
Phosphate  Bonded  Silicon  Nitride 
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As  discussed  in  Appendix  B,  a  series  solution  can  also  be  obtained  (instead  of  that  in 
terms  of  a  double  integral)  by  employing  the  calculation  of  residues.  However,  determina¬ 
tion  of  the  poles  of  the  integrand  (in  the  integral)  is  a  complicated  process,  when  several 
layers  are  involved.  The  poles  can  be  determined  easily  for  the  case  of  a  conductor  backed 
slab  of  one  material.  Some  results  of  computation  for  a  rectangular  wave  guide  are  given  in 
Appendix  C.  It  may  be  noted  that  with  an  accurate  integration  scheme  (some  times  up  to 
12,000  integration  points  are  required  for  one  integral),  the  results  from  the  two  approaches 
(integral  and  series)  are  identical  and  they  compare  reasonably  well  with  some  of  the  results 
reported  by  Bakhtiari  (Ref.  12).  However,  some  of  Bakhtiari’s  results  (Ref.  13)  do  not  agree 
with  those  obtained  using  the  series  and  the  integral.  It  is  likely  that  these  solutions  of  Bak¬ 
htiari  are  either  not  reported  correctly  or  the  error  may  be  the  result  of  inaccurate  evaluation 

of  the  integrals. 

Results  of  a  detailed  parametric  study  on  the  effects  of  air  gaps  (0  to  0.5  )  and  fre¬ 
quency  (8  to  12  GHz)  on  the  actual  conductance  and  susceptance  due  to  a  dielectric  mate- 
r-|al  _  2.8,  tan  8  =  0.05)  of  various  thicknesses  (0.075",  0.125",  0.235",  0.38  and  a 
half  space)  are  presented  in  Appendix  D.  The  results  show  that  for  slabs  of  small  thick¬ 
nesses  (0.075"  and  0.125"),  small  air  gaps  (up  to  0.1")  have  a  small  influence  on  the  con¬ 
ductance  for  a  frequency  of  8GHz,  but  the  susceptances  tend  to  decrease  rapidly  with  air 

gap  thickness. 

Dielectric  properties  of  layered  materials  are  often  influenced  by  air  gaps  or  delamina¬ 
tions,  which  may  be  present  between  the  layers.  Use  of  microwave  NDE  for  detecting  the 
air  gap  at  such  delaminations  has  been  suggested  by  some  investigators.  In  this  approach 
the  sensor  (a  rectangular  wave  guide)  is  placed  in  front  of  a  coupling  dielectric  sheet,  which 
is  backed  by  the  specimen  being  interrogated.  An  unintentional  air  gap  is  usually  present 
between  the  coupling  material  and  the  specimen  surface.  This  is  called  the  lift-off  effect 
and  one  first  attempts  to  find  the  direction  of  the  change  in  the  tip  of  the  vector  (complex 
reflection  coefficient)  due  to  small  lift-off  (air  gap).  This  direction  (that  of  the  change)  is 
called  the  lift-off  direction.  For  the  purpose  of  detecting  the  presence  or  absence  of  a  dela¬ 
mination  type  air  gap  inside  the  specimen,  the  change  in  the  component  of  the  complex  re¬ 
flection  coefficient  vector,  which  is  orthogonal  to  the  lift-off  direction,  is  normally  used,  so 
as  to  minimize  the  lift-off  effects  produced  by  wobbling  of  the  wave  guide  and  the  coupling 
sheet  over  the  specimen  during  scanning.  At  certain  frequencies  the  lift-off  is  a  minimum 
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and  the  orthogonal  component  is  a  maximum.  Appendix  E  describes  the  results  of  a 
parametric  study  to  find  the  frequency  ranges  which  should  be  used  so  as  to  minimize  the 
lift-off  effect. 

Appendix  F  gives  some  comparison  of  results  obtained  using  the  computer  code  ADMIT 
for  the  co-axial  line  with  the  results  presented  by  Bakhtiari  et  al  (Ref.  14)  and  Li  et  al  (Ref. 

1 5)  The  results  show  that  the  computer  code  is  free  from  any  bugs  and  the  results  can  be 
used  with  confidence  for  selecting  frequencies  and  other  parameters  for  interrogating 
specimens. 

In  Appendix  G,  an  attempt  is  made  to  correlate  experimental  observation  of  Liu  (Ref.  16) 
on  air  gap  and  other  effects  on  measurement  of  dielectric  properties  of  glass  fiber  reinforced 
composites  using  the  commercially  available  Hewlett  Packard  85070a  co-axial' probe  (Ref. 

1 7).  The  probe  appears  to  be  a  useful  tool  for  interrogating  specimens  for  detecting  varia¬ 
tions  in  dielectric  properties. 

TASK  4.  NDE  OF  SAMPLES 

Details  of  works  performed  at  VPI&SU  are  given  in  Part  2  of  the  report.  Phosphate 
Bonded  Silicon  Nitride  (PBSN)  disc  samples  prepared  for  the  program  (see  Task  1  described 
earlier  and  Part  2  of  the  report)  were  examined  using  ultrasonic,  infrared  (thermal),  optical 
reflection/transmission  and  microwave  imaging.  Transmitted  light  was  not  found  to  be  a 
viable  technique  for  inspection.  Internal  flaws  (such  as  voids)  could  be  detected  using  infra¬ 
red  technique  when  a  temperature  gradient  was  introduced  across  the  thickness.  Ultrasonic 
inspection  was  found  to  be  a  useful  tool  for  characterizing  flaws  due  to  porosity,  delamina¬ 
tion,  crack  and  metal  powder  inclusions  in  the  disc  samples.  Since  the  pores  reduce  the  ef¬ 
fective  modulus  of  the  material,  they  reduced  the  velocity  of  sound  in  PBSN  (Figure  20).  On 
the  other  hand,  when  immersed  in  water  for  about  four  hours,  the  pores  absorbed  water 
and  the  compressive  wave  speed  of  PBSN  increased  due  to  the  incompressibility  of  the  wa¬ 
ter  trapped  in  the  pores  (Figure  21).  Water  immersion  is  also  useful  since  water  acts  as  a 
good  coupling  agent  for  transmission  of  ultrasonic  waves.  Microwave  images  of  the  discs 
obtained  using  a  rectangular  wave  guide  (with  dimensions  0.9"  x  0.4"  and  a  frequency 
range  of  8.2  GHz  to  12.4  GHz)  are  shown  in  Figure  22.  It  may  be  noted  that  since  the 
wave  guide  dimensions  are  not  small,  there  is  distinct  edge  effect  near  the  free  edges  of  the 
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discs,  but  these  effects  should  not  be  of  importance  for  inspection  of  larger  samples  or  dome 
shaped  radomes. 


C-SCAN  PLOT  COLLECTED  ON:  Wed  Jan  24  1 7:24:34  1 996 


PART  I.D.  Number: 

disc  4 

PART  TYPE: 

disc 

Operator: 

Jon 

Transducer  Type: 

25  MHz 

Transducer  S/N: 

Transducer  Mfr: 

Panametrics 

Nominal  Cntr  Freq: 

25  MHz 

Attenuation: 

14 

Energy: 

4 

Damping: 

0 

Gain: 

60  dB 

HP  Filter: 

.03 

Commentl: 

1/24/95 

Comment2: 

Comment3: 

Comment4: 

Disc  was  constructed  to  have  a  porosity  gradient.  The  center  of  the  porosity 
"well"  appears  to  be  slightly  off  from  the  disc’s  center. 


Figure  20.  Ultrasonic  Image  of  a  Disk  #4  with  a  Porosity  Gradient 
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MICROWAVE  IMAGES  OF  DEFECTS  IN 
PHOSPHATE-BONDED  SILICON  NITRIDE 


TASK  5.  DATA  ACQUISITION  FROM  NDE  INSPECTION 


Based  on  the  evaluation  of  NDE  methods,  it  was  determined  that  a  reliable  assessment 
of  ceramic  composite  radomes  made  of  materials  similar  to  PBSN  would  require  ultrasonic 
pulse  echo  scanning  and/or  microwave  NDE.  A  commercially  available  scanning  package 
(hardware  and  software  from  Sonix)  capable  of  collecting  data  (ultrasonic  wave  form)  at 
digitization  rates  as  high  as  800  MHz,  using  a  variety  of  scan  procedures  while  tracking  the 
amplitude  and  time  of  flight  signals  within  eight  separate  gates,  was  selected.  The  hard¬ 
ware  was  adopted  to  a  custom  design  developed  at  VPI&SU  for  rotationally  symmetric  bul¬ 
let  (radome)  shaped  hollow  objects.  Another  software  package  capable  of  displaying  NDE 
data  at  a  moderate  level  of  resolution  (for  scans  at  moderately  high  speeds)  and  also  at  a 
higher  resolution  (for  more  careful  scanning  in  selected  areas)  was  utilized.  A  detailed  de¬ 
scription  of  data  collection  and  data  display  procedures  is  given  in  Sections  IV  and  V  of  Part 
2  of  the  report. 


c_r> 


Figure  23.  A  Schematic  Drawing  of  the  Inspection  System 
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TASK  8.  INSPECTION  PLAN 


Considerable  effort  was  directed  in  selecting  an  inspection  plan  that  would  combine 
economy,  flexibility,  and  reliable  application  of  ultrasonic  and/or  microwave  NDE  to  radomes 
of  different  sizes.  A  combination  of  commercially  available  computer  hardware  and  soft¬ 
ware,  fixtures  and  circuitry  units  designed  and  fabricated  at  VPI&SU,  as  well  as  various 
components  such  as  the  stainless  steel  holding  tank  available  from  NSWC,  were  utilized. 
Figure  23  gives  a  schematic  diagram  of  the  system.  The  NDE  unit  can  be  either  a  micro- 
wave  network  analyzer  or  an  ultrasonic  pulser/receiver.  For  ultrasonic  scanning,  the  dome 
supported  by  a  cradle  has  to  be  Immersed  in  water  in  a  holding  tank.  A  full  description  of 
the  holding  tank,  the  dome  cradle  and  rotation  system,  the  probe  holder  and  the  translation 
system,  the  motion  control  and  data  acquisition  system,  various  computing  hardwares  and 
softwares,  as  well  as  a  description  of  the  scanning  procedure  is  given  in  Part  2  of  the  re¬ 
port.  Some  images  which  were  also  obtained  for  a  full  scale  radome  supplied  by  NSWC. 
The  material  for  the  radome  was,  however,  not  PBSN.  Some  examples  of  these  images 
(peak  and  TOF)  with  inherent  impressions  and  seams,  etc.,  are  shown  in  Figures  24-27.  Im¬ 
ages  of  introduced  defects  (dents  and  cracks)  at  different  resolutions  are  illustrated  in  Fig¬ 
ures  28-31. 
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Figure  24.  Inherent  Impressions  and  Seams  (Complete  Image) 


Figure  25.  Impression  and  Seams  (Smaller  Area) 
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Figure  26.  Inherent  Impressions  (Peak) 


Figure  27.  Inherent  Impression  (TOF) 
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Figure  28.  Dent  (Peak) 


Figure  29.  Dent  (TOF) 
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Figure  30.  Crack  Tip  (Peak) 


Figure  31.  Crack  Tip  (TOF) 
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TASK  9.  SUGGESTIONS  FOR  ACCEPT/REJECT  CRITERIA 

The  results  in  Task  8,  summarized  earlier  and  described  in  detail  in  Part  2,  indicate  that 
the  inspection  system  is  capable  of  detecting  different  types  flaws  in  ceramic  composite  ra¬ 
domes.  In  particular,  a  map  of  the  radome  surface  showing  porosity  induced  variations  in 
the  ultrasonic  characteristics  (variations  in  sound  wave  velocity  or  time  of  flight)  can  be 
generated  using  the  inspection  system.  A  procedure,  described  in  Appendix  D  of  Part  2  of 
the  report,  allows  quantitative  assessment  of  the  level  of  porosity  for  the  Phosphate  Bonded 
Silicon  Nitride  material  used  in  this  study.  By  the  use  of  results  described  in  Task  2.1,  the 
variation  of  the  dielectric  constant  over  the  radome,  particularly  in  areas  or  patches  where  it 
is  maximum,  can  be  estimated.  Microwave  NDE  using  commercially  available  rectangular 
wave  guide  or  co-axial  cable  type  sensors  discussed  in  Task  3  should  also  be  performed  for 
direct  quantification  of  the  dielectric  constant  variation.  An  accept/reject  criterion  can  be 
chosen  based  on  some  acceptable  variation  (maximum  deviation  and  the  size  of  patch  de¬ 
fect)  of  porosity  and/or  dielectric  constant.  From  the  results  presented  in  Task  2.2,  it  ap¬ 
pears  that  a  dielectric  constant  variation  of  20%  over  a  single  patch  of  size  4"  x  4"  in  a 
typical  radome  (base  diameter  =  12.88",  height  =  26.13")  may  be  acceptable  provided  ap¬ 
propriate  corrections  are  made  to  account  for  expected  changes  in  boresight  errors  and  an¬ 
tenna  receiving  patterns.  Alternatively,  detailed  analytical  evaluation  of  the  effect  of  the 
patch  on  radome  performance  may  be  performed,  if  desired,  using  RACP  type  codes  used  in 
Task  2.2.  It  should  be  noted  that  correlations  of  analysis  results  with  some  experimental 
measurements  of  antenna  receiving  patterns  using  radomes  containing  patch  type  defects 
are  needed  before  such  results  are  used  in  practice. 

Since  materials  other  than  Phosphate  Bonded  Silicon  Nitride  are  also  being  considered 
for  next  generation  radomes  and  no  experimental  data  are  available  to  establish  a  quantita¬ 
tive  relationship  relating  the  porosity  to  the  ultrasonic  characteristics  and  also  the  dielectric 
constant,  some  relations  for  dilute  concentration  of  porosity  (or  its  change  in  a  material  with 
constant  porosity)  given  below  may  be  used  for  decision  making  purposes.  These  relations 
are  based  on  well  known  results  for  determining  the  effective  elastic  and  dielectric  con¬ 
stants  (Refs.  1 8,  9)  of  a  composite  made  of  an  isotropic  material  and  a  dilute  concentration 
of  spherical  voids  (or  inclusions).  Volume  fraction  of  voids  is  denoted  by  V  and  effective 
properties  are  denoted  with  an  asterisk. 
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Shear  Modulus 


v  being  the  Poisson's  ratio. 


Bulk  Modulus  (with  pores) 


Bulk  Modulus  (with  pores  filled  with 


’■A'-'- 

‘■-'['-'■sNs. 


water  treated  as  rigid  inclusions)  k*  =  k  1  +  V 


3(1 -y) 

(1+v) 


Density  with  pores  P  —  />  (l  —  l7) 

Density  with  pores  filled  with  water  p*  =  p  l-V  1- 

\  P  y 


Pw  is  the  density  of  water. 


The  wave  velocity  of  interest,  the  longitudinal  wave  velocity  in  the  composite,  is  given 


c;={(r+4///3)//>}'!  m 

Using  equations  (2)  -  (7)  and  neglecting  terms  with  higher  powers  of  V,  one  obtains  the  fol¬ 
lowing  expressions  for  changes  in  the  wave  velocity  due  dilute  concentration  of  pores.  For 
pores  not  filled  with  water 

c'-cfi  rf  1+v  I  5(1-2v)-1ll  (8) 

‘  T  K'-2  y)  7-5  v  2  J  J 

with  v  =  0.3 


Cl=Cj-0.676K] 

For  pores  saturated  with  water  after  long  immersion 

c*  =  c\ 1 + v\i  -  -  --^4 

1  [  7  -5v  2  p  Jj 


with  v  =  0.3  and  pw/p  =  0.4 
Cl=CLl 1-0.436T] 


(ID 


Equations  (9)  and  (11)  indicate  that  for  materials  with  Poisson's  ratio  and  density  compara¬ 
ble  to  those  listed  before  these  equations,  the  wave  velocity  decreases  by  3.4%  due  to  5% 
void  content,  whereas  when  saturated  with  water  it  increases  by  2.2%. 

It  may  be  noted  that  the  test  data  for  unsaturated  PBSN  samples  with  voids,  given  in 
Appendix  D  of  Part  2,  do  not  correlate  with  equation  (9),  i.e.,  the  change  in  velocity  appears 
to  be  much  higher  than  that  given  by  (9).  There  are  two  possible  reasons  for  this  discrep¬ 
ancy,  namely,  (i)  the  pores  may  not  be  spherical  (it  is  known  (Ref.  1 9)  that  for  oblate  sphe¬ 
roidal  pores  which  are  crack  like,  the  effective  moduli  decrease  at  a  faster  rate  with  V)  and 
(ii)  the  measured  porosity  change  is  that  of  connected  porosity,  which  may  be  lower  than 
the  total  change.  Moreover,  the  pores  may  contain  some  solids,  which  may  also  reduce  the 
change  in  the  effective  densities  given  by  equations  (5)  and  (6).  It  appears,  therefore,  that 
a  better  estimate  may  be  obtained  by  neglecting  the  density  change  and  multiplying  the 
change  in  the  velocity  by  a  factor  a,  which  yield  the  following  estimates.  For  pores  not 
filled  with  water 


c;=cL 


1  -aV\ 


1+v 
4(l -2v) 


5(i-2k)] 
7-5v  J 


(8a) 


For  v  =  0.3 

C;=Q[l-1.18av] 

For  pores  saturated  with  water 


r*  -C 


1- 0.5- 


5(1 -2^)1 

7-5v  J 


(9a) 


(10a) 


For  v  =  0.3 

C*L  =  CL[l  +  0. \4aV]  (11a) 

From  the  test  data  given  in  Part  2,  it  appears  that  a  in  (9a)  must  be  of  the  order  of  4  for 
PBSN.  It  may  be  noted  that  for  oblate  (flat)  spheroidal  pores  with  aspect  ratio  of  the  order 
of  0.1,  the  decrease  in  effective  moduli  may  be  three  to  four  times  that  for  spherical  pores 
(Ref.  19). 

Equations  (8a)  and  (10a)  with  a«  4  can  be  used  to  estimate  the  small  variations  of  po¬ 
rosity  (much  larger  variations  will  possibly  be  not  acceptable)  from  ultrasonic  characteriza¬ 
tion,  although  it  may  be  useful  to  perform  some  tests  for  new  material  systems  (with  known 
porosities)  under  development,  in  order  to  obtain  a  correct  estimate  of  a. 
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The  changes  in  the  dielectric  constant  due  to  small  variation  of  porosity  can  be  esti¬ 
mated  from  the  following  equations 


(12) 

where  ea  =  1.0,  the  dielectric  constant  of  air.  However,  since  the  dielectric  solid  is  some¬ 
times  designed  to  have  a  constant  porosity  which  may  be  high  (as  in  PBSN)  and  the  dielec¬ 
tric  constant  for  the  solid  itself  may  be  quite  high,  a  more  accurate  expression  (Ref.  9)  is 
given  by 

e  =£s 

where 

es  =  the  dielectric  constant  for  the  solid  without  any  voids 
V,  =  total  void  content  inclusive  of  desired  one 


1- 


K-3  kzld 

ea+2es+Vt(es-eJi 


(13) 
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V-3(s-e.) 


2  e  +  e. 


a  J 


Equation  (13)  is  an  upperbound  result  for  nondilute  concentration  of  spherical  voids  and  has 
been  used  in  Task  2.1  to  estimate  the  effects  of  porosity  on  the  dielectric  constant.  Effects 
of  variations  in  porosity  on  the  change  in  the  loss  tangent  of  the  material  can  be  neglected 
since  the  loss  is  usually  small.  No  estimates  are  currently  available  for  the  effect  of  the  as¬ 
pect  ratio  of  oblate  spheroidal  voids  on  the  effective  dielectric  constant.  It  is  expected, 
however,  if  s  (or  8S)  in  (12)  (or  (13))  are  not  too  large  these  relations  should  yield  a  good  es¬ 
timate  of  the  change  in  e\ 

Once  the  variations  in  the  dielectric  constant  over  a  radome  of  next  generation  material 
are  estimated  from  ultrasonic  characterization,  it  would  be  worthwhile  to  use  microwave 
NDE  for  obtaining  a  direct  quantification  and  correlation.  No  detailed  microwave  evaluations 
of  the  PBSN  material  were  performed  in  this  program,  but  commercially  available  probes, 
such  as  the  Hewlett  Packard  85070A  and  the  8722C  vector  network  analyzer  may  be  util¬ 
ized,  as  discussed  in  Task  3  and  Ref.  16.  Based  on  the  estimated  changes  in  the  dielectric 
constant,  a  decision  regarding  the  acceptance  of  a  radome  can  be  made,  if  some  limits  are 
chosen  in  the  dielectric  property  variation  as  described  earlier.  Analytical  evaluation  of  the 
electromagnetic  performance  of  radome  can  also  be  made  to  determine  the  acceptability  of 
a  radome. 
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Effects  of  porosity  and/or  surface  crack  type  defects  on  stiffness,  strength  and  thermal 
expansion  characteristics  can  be  estimated  using  the  methods  and  results  given  in  Task  2.1. 
However,  dilute  concentration  of  additional  porosity  over  the  desired  value  and  small  crack 
like  defects  (sizes  less  than  those  described  in  Task  2.1)  are  not  expected  to  have  any  sig¬ 
nificant  influence  on  thermomechanical  properties  of  the  material. 
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CONCLUDING  REMARKS 


Optical  reflection/transmission  measurements  used  for  nondestructive  evaluation  of  cur¬ 
rently  used  radomes  are  not  suitable  for  use  with  next  generation  ceramic  composite  ra¬ 
domes,  since  these  materials  do  not  allow  light  transmission.  Ultrasonic  characterization 
has  been  found  to  be  a  useful  tool  for  quantification  of  porosity  variations  and  other  de¬ 
fects,  such  as  crack  like  defects.  Variation  in  porosity  and/or  reinforcements  is  the  domi¬ 
nant  anomaly  in  ceramic  matrix  radomes  and  it  influences  the  dielectric  property  of  the  ma¬ 
terial,  which  in  turn  changes  the  electromagnetic  performance  of  the  radome.  An  inspection 
system  for  interrogation  of  radome  shaped  shells  has  been  developed.  The  system  can  ac¬ 
commodate  radomes  of  various  sizes  and  has  been  shown  to  generate  a  map  of  property 
variations  (especially  the  change  in  porosity)  using  ultrasonic  probes.  Methods  of  estimating 
porosity  variation  from  ultrasonic  characteristics  and  the  changes  in  the  dielectric  constant 
from  those  in  porosity,  which  can  be  used  for  next  generation  ceramic  composites,  are  de¬ 
scribed.  Approaches  for  selecting  accept/reject  criteria  based  on  the  results  of  inspection 
are  suggested,  which  should  be  examined  in  conjunction  with  some  trial  experiments  and 
correlation  of  test  results  with  analytical  evaluation  of  the  electromagnetic  performance  of 
radomes  made  with  intentional  porosity  variations  over  a  patch.  The  RACP  code  has  been 
modified  for  analytical  evaluation. 

The  inspection  system  can  also  make  use  of  microwave  probes,  although  a  demonstra¬ 
tion  of  such  imaging  could  not  be  performed  on  full  scale  radomes.  Microwave  probes  can 
yield  a  direct  measure  of  dielectric  property  changes  and  was  used  for  some  small  samples. 
However,  edge  effects  influence  the  image  for  small  samples  since  the  probes  are  usually  of 
moderate  sizes.  Such  edge  effects  are  not  expected  in  imaging  of  full  scale  radomes.  There¬ 
fore,  microwave  imaging  of  radomes  using  the  inspection  system  should  be  attempted  in  the 
future.  However,  microwave  probes  need  a  conducting  flange  and  a  customized  design  of 
the  flange  may  be  required  for  the  curvature  in  radomes.  Attempts  were  made  for  analytical 
evaluation  of  unavoidable  air  gaps  between  the  flange  and  the  sample  being  interrogated 
(due  to  surface  roughness  and/or  curvature)  as  well  as  other  parameters  (such  as  interroga¬ 
tion  frequency)  on  expected  signals  and  the  results  are  discussed. 

Attempts  are  being  made  to  contact  some  manufacturers  of  next  generation  radomes 
and  to  determine  their  interest  in  the  inspection  system  developed  in  the  program. 
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INTRODUCTION 


The  RACP™  code  is  quite  versatile  and  allows  evaluation  of  boresight  errors  for 
different  types  of  radomes  with  variable  thickness  and  properties  as  well  as  monopulse 
antenna  near  fields  (receiving  patterns)  for  various  polarizations  of  the  incident  plane  wave 
and  the  antenna.  However,  in  its  present  form  the  code  does  not  allow  consideration  of  a 
patch  type  defect  in  the  radome  (Figure  A1).  The  code  was,  therefore,  modified  to 
consider  such  property  variations  over  a  region  bounded  by  two  prescribed  axial 
coordinates  (z)  and  circumferential  coordinates  (<|>)  as  shown  in  Figure  A2.  Since  it  is 
expected  that  the  property  change  will  be  maximum  at  the  center  of  the  patch  and  will 
reduce  to  zero  at  its  boundaries,  it  is  assumed  that  the  property  change  (real  and  imaginary 


part  of  the  dielectric  constant  or  thickness)  is  of  the  following  form 

.  .  7t(z  -  Z\) 

A  =  Ao  cos— - -cos  ^  - 

e  e 


(A-1) 


where  z\  and  <j>i  are  the  coordinates  of  the  center  of  the  patch  and  «  and  6  quantify  its 
dimension  in  axial  and  circumferential  direction.  It  may  be  noted  that  <|>i  defines  the 
relative  angular  position  of  the  patch  with  respect  to  the  antenna.  <|>i  =  90°  when  an  axial 
plane  (through  radome  z-axis)  passing  through  the  center  of  the  . patch  coincides  with  the 
azimuth  plane  (xAzA)  of  the  antenna  (Figure  A2).  For  <|>,  =  0°,  the  axial  plane  coincides  with 
yAzA  or  elevation  plane  of  the  antenna  (Figure  A2). 

An  illustrative  example  and  some  results  are  given  in  pages  A3-A18. 

The  changes  in  input  data  for  RACP™  required  to  analyze  the  effects  of  patch  defects 
are  described  in  pages  A19-A24. 
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Axial  Distance  (In.) 


Figure  A1 .  Radome  with  a  Patch  Defect,  p.  8x 


Figure  A2.  Dielectric  Property  Validation  and  Position  of  the  Patch,  p.8y 


AN  ILLUSTRATIVE  EXAMPLE 


A  typical  tangent  ogive  radome  shown  with  the  antenna  in  Figure  A1  was  considered 
for  demonstration  purposes.  The  following  geometric  and  material  properties  were  used 
(Figure  A3). 

Radome  base  diameter,  Dos  =  1 2.88  inch 

Radome  length,  Los  =  26.43  inch 

Bulkhead  z  coordinate,  zBot  =  0.0  inch 

Antenna  Aperture  Radius,  R»p  =  5.0  inch 

Gimbal  point  from  radome  base,  Rr  =  7.057  inch 

Antenna  aperture  offset  from  gimbal,  Ra  =  1 .908  inch 

Frequency,  F  =  1 1 .803  GHZ 

Radome  wall  thickness,  DIN  =  0.351  inch 

Dielectric  constant  er,  ER  =  5.0 

Loss  tangent,  TD  =  0.03 

For  estimation  of  changes  in  boresight  errors  due  to  the  patch,  the  dielectric  constant,  e,, 
at  center  was  assumed  to  be  4.0  (Figure  A2).  For  estimation  of  changes  in  receiving 
patterns  one  more  value  of  s,  at  the  center  of  the  patch  (sr  =  2.0)  was  also  selected. 

Variations  of  boresight  errors,  error  slopes  and  relative  antenna  gain  with  yaw  gimbal 
angle  for  different  patch  sizes  (at  an  axial  position  z,  =  15"  and  circumferential  position  <J>i 
=  90°),  when  the  antenna  and  the  incident  wave  are  both  vertically  polarized,  are  shown 
in  Figures  A4  -  A8.  Variations  of  the  same  quantities  with  pitch  gimbal  angle  are  plotted  in 
Figures  A9  -  A13.  It  may  be  noted  that  the  elevation  errors  (and  slopes)  are  zero  (Figures 
A4  and  A6)  for  this  case.  Also,  the  azimuth  error  increases  with  increasing  patch  size  for 
positive  yaw  gimbal  angles  (zA  of  antenna  moves  towards  the  patch). 

Effects  of  angular  position  (<|>  =0  to  90°)  for  a  6"  patch,  when  the  antenna  is  right 
hand  circularly  polarized  and  the  incident  wave  is  vertically  polarized,  are  shown  in  Figures 
A14  -  A23.  Boresight  elevation  errors  (and  slopes)  are  not  zero  (Figures  A14  and  A16)  in 
this  case. 


A-3 


Computed  amplitude  and  phase  of  the  antenna  receiving  patterns  at  elevation  and 
azimuth  cuts  for  a  6"  patch  defect  (with  s r  =  4.0  and  2.0  at  the  center)  located  at  axial 
position  z\  —  15"  and  circumferential  position  4>i  =  90°  for  one  value  of  the  yaw  gimbal 
angle  (=  10°)  are  shown  in  Figures  A24  -  A31 . 


Figure  A3.  Geometry  and  Dimensions  of  Radome  and  Antenna 
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Bore  S)gM  A**nulh  Error  Stop*,  (%> 


.20  -15  -10  -6  0  6  10  IS  20 

Yew  Glmbel  Angle  There,  (degree*) 

Figure  A6.  Effect  of  Defect  Size  at  90  Degree  Position  and  z  =  1 5  in.  on  Bore  Sight 
Elevation  Error  Slope  for  Vertically  Polarized  Antenna  and  Incident  Wave 


Yew  Gmtfcnf  Angle  Theta,  (degrees) 


Figure  A7.  Effect  of  Defect  Size  at  90  Degree  Position  and  z=  15  in.  on  Bore  Sight 
Azimuth  Error  Slope  for  Vertically  Polarized  Antenna  and  Incident  Wave 
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Yaw  Gfanbaf  Angle  TheU,  (Agrees) 


Figure  A8.  Effect  of  Defect  Size  at  90  Degree  Position  and  z  =  1 5  in.  on  Relative  Antenna 
Gain  for  Vertically  Polarized  Antenna  and  Incident  Wave 


Pitch  Grmbal  Angle  The**,  {degrees) 


Figure  A9.  Effect  of  Defect  Size  at  90  Degree  Position  and  z=  15  in.  on  Bore  Sight 
Elevation  Error  for  Vertically  Polarized  Antenna  and  Incident  Wave 
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Bor*  Sight  Elevation  Error  Slop*,  <%) 
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Figure  A 10.  Effect  of  Defect  Size  at  90  Degree  Position  and  z  =  1 5  in.  on  Bore  Sight 
Azimuth  Error  for  Vertically  Polarized  Antenna  and  Incident  Wave 
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■2  Inch  Defect  Size 
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Figure  A1 1 .  Effect  of  Defect  Size  at  90  Degree  Position  and  z  =  1 5  in.  on  Bore  Sight 
Elevation  Error  Slope  for  Vertically  Polarized  Antenna  and  Incident  Wave 
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Figure  A1 2.  Effect  of  Defect  Size  at  90  Degree  Position  and  z  =  1 5  in.  on  Bore  Sight 
Azimuth  Error  Slope  for  Vertically  Polarized  Antenna  and  Incident  Wave 


■  ♦—<)  heft  Defect  Size  : 
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Figure  A1 3.  Effect  of  Defect  Size  at  90  Degree  Position  and  z  =  1 5  in.  on  Relative  Antenna 
Gain  for  Vertically  Polarized  Antenna  and  Incident  Wave 
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Figure  A14.  Effect  of  Angular  Position  of  6  in.  Defect  at  z  =  1 5  in.  on  Bore  Sight  Elevation 
Error  for  RHC  Polarized  Antenna  and  Vertically  Polarized  Incident  Wave 
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—90.0  degree.  6  inch  Defect 
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Figure  A1 5.  Effect  of  Angular  Position  of  6  in.  Defect  at  z=  1 5  in.  on  Bore  Sight  Azimuth 
Error  for  RHC  Polarized  Antenna  and  Vertically  Polarized  Incident  Wave 
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Figure  A1 6.  Effect  of  Angular  Position  of  6  in.  Defect  at  z  =  1 5  in.  on  Bore  Sight  Elevation 
Error  Slope  for  RHC  Polarized  Antenna  and  Vertically  Polarized  Incident  Wave 


Yaw  Gknt>al  Angle  Theta,  (degrees) 


Figure  A1 7.  Effect  of  Angular  Position  of  6  in.  Defect  at  z  =  1 5  in.  on  Bore  Sight  Azimuth 
Error  Slope  for  RHC  Polarized  Antenna  and  Vertically  Polarized  Incident  Wave 
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Figure  A18.  Effect  of  Angular  Position  of  6  in.  Defect  at  z=15  in.  on  Relative  Antenna 
Gain  for  RHC  Polarized  Antenna  and  Vertically  Polarized  Incident  Wave 


— No  Oefect 
— 90.0  degree,  6  Inch  Defect 
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Figure  A19.  Effect  of  Angular  Position  of  6  in.  Defect  at  z  =  1 5  in.  on  Bore  Sight  Elevation 
Error  for  RHC  Polarized  Antenna  and  Vertically  Polarized  Incident  Wave 
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Figure  A20.  Effect  of  Angular  Position  of  6  in.  Defect  at  z  =  1 5  in.  on  Bore  Sight  Azimuth 
Error  for  RHC  Polarized  Antenna  and  Vertically  Polarized  Incident  Wave 
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Figure  A21 .  Effect  of  Angular  Position  of  6  in.  Defect  at  z=  1 5  in.  on  Bore  Sight  Elevation 
Error  Slope  for  RHC  Polarized  Antenna  and  Vertically  Polarized  Incident  Wave 
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Figure  A22.  Effect  of  Angular  Position  of  6  in.  Defect  at  z  =  1 5  in.  on  Bore  Sight  Azimuth 
Error  Slope  for  RHC  Polarized  Antenna  and  Vertically  Polarized  Incident  Wave 
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Figure  A23.  Effect  of  Angular  Position  of  6  in.  Defect  at  z  =  1 5  in.  on  Relative  Antenna 
Gain  for  RHC  Polarized  Antenna  and  Vertically  Polarized  Incident  Wave 
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PHASE*  DECREE.  DELTA  ELEVATION  ELEVATION  COMPONENT 


Figure  A26.  Effect  of  Defect  on  Receiving  Pattern-Delta  Amplitude  Elevation  Cut 
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Figure  A27.  Effect  of  Defect  on  Receiving  Pattern-Delta  Phase  Elevation  Cut 
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Figure  A28.  Effect  of  Defect  on  Receiving  Pattern-Sum  Channel  Amplitude  Azimuth  Cut 
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Figure  A29.  Effect  of  Defect  on  Receiving  Pattern-Sum  Channel  Phase  Azimuth  Cut 
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Figure  A30.  Effect  of  Defect  on  Receiving  Pattern-Delta  Amplitude  Azimuth  Cut 
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Figure  A31 .  Effect  of  Defect  on  Receiving  Pattern-Delta  Azimuth  Amplitude  Azimuth  Cut 
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CHANGES  IN  INPUT  REQUIRED 


The  input  data  file  for  RACP™  is  shown  in  Table  A1 .  The  case  of  patch  type  defect  is 
chosen  through  the  advanced  data  feature  by  setting  SETDEF  =  F  and  the  logical  variables 
LVTHK  and/or  LERTD  eoual  to  T  (True)  and  the  following  additional  inputs  to  describe  the 
patch  geometry  and  the  variations  in  layer  thicknesses  and/or  dielectric  properties  as 
described  below.  Note  that  LVTHK  and  LERTD  are  used  in  the  present  version  of  RACP™ 
for  prescribing  changes  in  layer  thicknesses  and  dielectric  properties,  which  vary  following 
a  power  law  in  axial  <z  direction)  and  as  a  sine  function  (Abs  (sin  (n<j>/2+\|/d)))  in  <}>  direction. 
The  value  n  (input  as  ND  or  NERTD)  prescribes  the  number  of  half  cycles  of  the 
circumferential  variation.  To  differentiate  between  such  a  variation  from  that  of  a  patch 
defect  ND  and/or  NERTD  are  set  equal  to  an  integer  greater  than  99  (  =  100  say). 

To  consider  a  thickness  variation  over  a  patch  described  by  equation  (A-1),  the 
following  inputs  are  required. 


LVTHK 

TZIN 

ND 

DDZIN(I) 


ZEXP(I) 


T 

Real  variable  equal  to  the  axial  radome  coordinate  of  the  center  of  the 
patch  in  inches,  i.e.,  the  value  of  z\  in  equation  (A-1 ). 

Integer  set  >  99  identifying  patch  defect. 

Real  array  (single  value  for  one  layer  thick  wall)  containing  the  change  in 
thickness  (in  inches)  of  layer  I  at  the  center  of  the  patch  from  the 
prescribed  thickness  DIN(I). 

The  array  of  ratios  z,/C(I),  e(I)  being  the  size  of  the  patch  in  inches  at 


layer  I 

(see  equation  (A-1 ))  in  axial  direction. 

DDCIN(I)  -  The  array  with  circumferential  dimensions  0  (equation  (A-1))  in  radians 
of  the  patch  in  layer  I. 

PSID(I) -  The  array  with  the  values  of  <j>i  in  radians  for  layer  I,  denoting  the 


position  of  the  center  of  the  patch. 


When  thickness  variations  are  absent  in  the  patch  LVTHK  may  be  set  equal  to  F  and  other 
inputs  are  not  required.  Alternatively  DDZIN  may  be  set  equal  to  0.0. 
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For  variations  in  dielectric  properties  over  a  patch  of  the  form  given  by  equation  (A-1), 
the  required  inputs  are  as  described  next. 


T 

Real  variable  equal  to  the  axial  radome  coordinate  of  the  center  of  the 
patch  in  inches,  i.e.,  the  value  of  zi  in  (A-1). 

Integer  set  >  99  identifying  patch  defect. 

Real  array  (single  value  for  one  layer)  containing  the  change  in  the 
dielectric  constant  of  layer  I  at  the  center  of  the  patch  from  the 
prescribed  value  ER  (A-1 ). 

The  array  of  ratios  zi/t(I),  {(I)  being  the  size  of  the  patch  in  inches  of 
layer  I  in  axial  direction  (see  (A-1)). 

The  array  with  circumferential  dimensions  0  (equation  (A-1 ))  in  radians 
of  the  patch  in  layer  I. 

The  array  with  the  values  of  4>i  in  radians  for  layer  I,  denoting  the 
position  of  the  center  of  the  patch. 

TDZ(I),  ZEXPTD(I),  TDC(I),  and  PSITD  (I)  -  These  quantities  describe  the  variations 
of  loss  tangents  in  the  same  way  as  ERZ(I),  ZEXPER(I),  ERC(I)  and 
PSIDER(I)  describe  the  variation  of  the  dielectric  constants  described 
above. 

NOTE:  It  is  important  to  note  that  it  is  advantageous  to  use  a  raster  scan  system  (by  using 
either  RASTER.COM  (with  geometric  optics  approximation)  or  RASTERNG.COM  (for  no 
geometric  optics  approximation)  as  the  CMD.COM  file  for  consideration  of  patch  type 
defects.  This  is  done  by  setting  COPY  RASTER.COM  (or  RASTERNG.COM)  CMD.COM 
before  execution.  Descriptions  of  the  raster  scan  system  and  the  polar  scan  system  are 
given  in  Chapter  4  of  the  RACP™  reference  library  (Ref.  20  in  the  text).  The  raster  scan  is 
commonly  used  in  RACP™  when  the  antenna  is  offset  from  the  axis  of  the  radome  defined 
by  the  coordinates  (Xo,  yc,  Zo  shown  in  Figure  A32  or  XROSP,  YROSP  and  ZROSP),  but  they 
may  be  set  to  zero. for  no  offsets.  The  angles  <J>P  and  0i  (Figure  A-32)  specify  the  plane  of 
scan  and  the  angle  (gimbal  angle)  in  the  scan  plane  of  the  radome  coordinate  system.  <)>„  = 
90°  specifies  the  elevation  yAzA  plane  (pitch  plane)  of  the  antenna.  0L  =  0  aligns  the  zR  (of 


LERTD  - 
TKIN 

NERTD  - 
ERZ(I) 

ZEXPER(I)  - 

ERC(I) 

PSIDER(I)  - 
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radome)  with  zA  (of  antenna).  Positive  0L  scans  the  radome  axis  in  positive  yA  direction.  <t»p 
=  constant  selects  a  great  circle  plane  that  is  approximately  vertical.  A  special  provision 
permits  a  single  horizontal  scan  (in  xAzA  plane  of  antenna,  i.e.,  yaw  or  azimuth  plane)  by 
setting  <j>p  =  0.0  (see  page  40  of  Ref.  20  of  text).  When  raster  scan  is  used,  the  angles 
PSID(I),  PSIDER(I),  and  PSITD(I)  defining  the  patch  circumferential  position  can  be  kept 
fixed  for  any  value  of  ^  or  0L.  On  the  other  hand,  if  polar  scan  is  used,  these  quantities 
have  to  be  changed  depending  on  the  value  of  <J>P,  the  plane  of  scan. 


Figure  A32.  Polar  Scan  and  Raster  Scan  Coordinate  Systems  Used  by  Program  RACP 
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I  f 


A  sample  input  data  set  for  a  6"  patch  defect  t  —  6"  with  center  located  at  z\  —  15 
(Zi/t  =  1 5/6  =  2.5)  and  4>i  =  0°  are  shown  in  Table  A2.  Note  that  the  radius  of  the 
radome  at  zi  =  15"  is  close  to  4"  so  that  the  value  of  0  (the  circumferential  dimension  of 
the  patch  is  equal  to  6/4)  is  1 .5  radians.  The  inputs  described  earlier  are  underlined  in 
Tables  A1  and  A2.  Changes  in  thickness,  dielectric  constant  and  loss  tangent  of  a  single 
layer  are  set  equal  to  zero  so  that  the  results  would  be  identical  to  the  case  with  no 
defects. 

Changes  have  been  made  in  the  subroutine  RXMIT  to  consider  the  effects  of  the  patch 
defect  by  the  use  of  relations  of  the  type  given  by  equation  (A-1)  for  computation  of  the 
thickness,  the  dielectric  constant  and/or  the  loss  tangent,  when  the  point  under 
consideration  falls  inside  the  patch.  It  should  be  noted  that  the  subroutines  SHAPE  or  SHP 
computes  the  radome  wall  thicknesses  for  plotting  purposes  only.  No  changes  have  been 
made  in  these  subroutines  so  the  shape  plot  obtained  may  be  erroneous  if  thickness 
variations  are  used  (LVTHK  set  equal  to  T),  but  boresight  errors  and  receiving  patterns  will 
be  correct. 
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Table  A1 .  Input  Data  File  for  RACP 


READ (5,6)  TITLE  (1X,A80  FORMAT) 

READ (5/ * )  SETDEF , GRAFS A, GRAFNF , GRAFRV , SUPPRS , IPENCD 

READ(5, *)  NFINE , NPHI , NTHE 

READ ( 5 , * )  DIAOS,RR, RA, RAPMAX 

IF  ( IPENCD .LE . 1 )  READ ( 5 , * )  FREQ,OSANG 

IF  (IPENCD. EO. 2)  READ ( 5 , * )  FSTART , FSTOP , DELF , OSANG 


READ (5, 


r) 


IOPT , I POL , ICASE , N , IPWR 
READ(5,*)  DIN(I) ,ER(I) ,TD(I) 

READ ( 5 , * )  FINR( I ) , ZTOPIN ( I ) 

READ ( 5 , * )  PHI ( I ) 

READ (  5  ,  *  )  THETA(I) 

READ ( 5 / * )  NX,NY,NREC 

READ (5,*)  NXE,MX, IMOD 

READ ( 5 , * )  ANG1, ANG2 

READ (5,*)  GRAFXN , GRAFSH ,  NZ 

READ ( 5 , * )  APIN, ZBOTIN, ISHAPE,ALFA 

READ (5,*)  IMTR,KMTR 

Read  the  following  advanced  feature  data  if  SETDEF=F: 


(I*=l  to  N  ) 

(1=1  to  NFINE) 

(1=1  to  NPHI  planes) 
(1=1  to  NTHETA) 


READ (5,1)  DUMMY 

READ ( 5  r  * )  XROSP , YROSP  , ZROSP 

READ (5, *,ERR=350)  LVTHK , LERTD 

IF  (LVTHK)  THEN 

READ ( 5 , * )  TZIN,ND 

READ  (  5  ,  * )  DDZINm,ZEXPm  ,DDCIN  ( I )  ,PSID  ( I )  (1=1, N) 

ENDIF 

IF  (LERTD)  THEN 

READ ( 5  ,  *  )  TKIN  ,  NERTD 

READ ( 5  ,  * )  ERZ ( I }  , ZEXPER (1) , ERC (II# PS IERD ( I )  (1=1, N) 

READ ( 5 , * )  TDZ (1) , ZEXPTD (II, TDC ( I ) ,PSITDD(I)  (1=1, N) 
ENDIF 

READ ( 5 , * , ERR=3 5 0 )  LATHK , LAERTD , GRAFAT , LPRTAT , LRDXMN 
IF  (LATHK)  THEN 

READ { 5 ,  * )  NPINP 

READ (5,*)  ZD ( I ) / ( DARB ( I ,  J ) , J=l,N)  (1=1, NPINP) 


ENDIF 


350  CONTINUE 


IF  ( ( . NOT . LVTHK ) .AND. ( .NOT. LERTD) .AND. ( .NOT. LATHK) .AND. 
$  ( . NOT . LAERTD ) )  THEN 
READ (  5  ,  *  , ERR=355 )  LRDXMN 
355  CONTINUE 
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Table  A2.  Sample  Input  Data  Set  for  6"  Defect 


RCAS  1  FOR  RACP  I SHAPE =3 

(SINGLE  LAYER 

TEST  CASE) 

F  F  F  F  T 

0 

1 

2  11 

12.8800000 

7.0570000 

1.9080000 

5.0000000 

11.8030000 

2.0000000 

1 

1  1 

1 

0 

0-700000E-001  S. 0000000  3.000000E-003 

2.0520190  26.4300000 

.0000000 
90.0000000 
.0000000 
2.0000000 
4.0000000 
6.0000000 
8.0000000 
10.0000000 
12.0000000 
14.0000000 
16.0000000 
18.0000000 
20.0000000 

16  16  32 

256  16  4 

-15.  15. 

T  T  50 

0.  0.  3  0. 

50  50 

THIS  LINE  AND  SUBSEQUENT  DATA  ARE  READ  IF  SETDEF=F . 

0.  0.  0. 

T _ T 

15. _ 100 

0.0 _ 2_;_5 _ 1.50  0. 


15. 

100 

0.0 

2.5 

1.50 

0. 

0.0 

2.5 

1.50 

0. 

F 

F  F  F 

F 
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Advanced 

Composites 

Technology 


introduction 


Flanged  open-ended  rectangular  wave  guides  and  co-axial  cables  are  being  employed  by 
many  investigators  for  microwave  NDE  of  ceramics  and/or  other  materials  used  in  various 
applications.  Such  interrogations  are  especially  suitable  for  cases  where  the  electro¬ 
magnetic  properties  of  the  materials  or  the  components  have  to  meet  certain  specific  re¬ 
quirements  (for  example,  in  radome  applications).  Dielectric  properties  of  ceramics  and  rein¬ 
forced  ceramics  are  often  strongly  influenced  by  the  reinforcement  and  void  content  as  well 
as  delaminations  and  other  types  of  flaws.  Interrogation  of  such  materials  to  determine 
such  flaws  and  anomalies  by  using  open-ended  sensors  appears  to  be  an  attractive  NDE 
technique,  where  a  probe  is  inserted  through  a  slot  or  other  opening  at  the  center  of  the 
wave  guide  (or  co-axial  cable)  to  monitor  the  variations  of  the  standing  wave.  Since  various 
parameters  other  than  the  dielectric  properties  of  the  materials  being  interrogated  (such  as 
air  gaps,  metallic  backing  plates  or  other  backing  materials  as  well  as  delaminations)  affect 
the  standing  wave  pattern,  analytical  estimates  of  the  expected  response  are  very  useful  for 
planning  the  experiments  and  for  correct  interpretation  of  the  measurements.  Various 
analyses  have  been  performed  for  this  purpose  (see,  for  example,  ref.  1 ,  2).  In  what  follows 
the  problem  shown  in  Figure  B-1  is  reconsidered  for  the  purpose  of  this  study. 


Metal  Backing 
Or  8acked  8y 
Material  Medium  (n  + 1 ) 


Figure  B-1 .  Layered  Medium 
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STATEMENT  OF  THE  PROBLEM 


For  interrogation,  the  dominant  TEM  mode  signal  from  the  wave  guide  hits  the  aperture 
area  over  the  layered  medium.  For  a  rectangular  wave  guide,  the  aperture  is  rectangular, 
whereas  for  a  co-axial  cable,  it  is  bounded  by  two  concentric  circles.  Outside  the  aperture, 
the  electric  fields  tangential  to  the  plane  of  the  aperture  (z  =  0,  Figure  B-1)  may  be  as¬ 
sumed  to  be  zero  because  of  the  presence  of  the  conducting  flange.  The  electric  and  mag¬ 
netic  vector  fields,  E  and  H  (a  bar  under  a  quantity  indicates  a  vector),  in  each  of  the  layers 
can  be  expressed  in  terms  of  the  tangential  components  of  the  electric  or  magnetic  fields  at 
the  aperture.  To  do  this,  the  solution  of  Maxwell's  equation  for  each  layer  is  sought  in 


terms  of  the  vector  potential  F  satisfying 

V2F  +  k?F  =  0  (B-1) 

F=Fxax+Fyay  (B-2) 

ax  ay  are  the  unit  vectors  in  x  and  y  directions.  E  and  H  are  given  by 

E  =  -VxF  (B-3) 

H  =— [k?F  +  VVF]  (B-4) 

\W 

where  j  =  V- 7  and  0)  is  the  frequency,  jx  is  the  magnetic  permeability  and  ki  is  the  com¬ 
plex  propagation  constant  of  the  layer  material,  i.e., 

kf-kj«-k2  (*’-■«')  (B-5) 


s',  s"  are  the  real  and  imaginary  parts  of  the  dielectric  constant  s  and  k<,  is  the  free  space 
propagation  constant  given  by 


A0  and  c  are  the  free  space  wave  length  and  wave  speed,  respectively. 
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SOLUTION  IN  TERMS  FOURIER  TRANSFORMS 


All  quantities  appearing  in  equations  (B1-B4)  can  be  assumed  as  a  product  of  two  func¬ 
tions,  one  varying  with  space  (x,  y,  z  coordinates)  and  the  other  with  time,  t.  The  time 
varying  function  is  the  same  for  all  quantities  and  is  assumed  to  be  exp  ( jcot ).  In  the  ex¬ 
pressions  of  the  quantities  which  follow,  this  time  dependent  function  will  be  omitted,  since 
their  spatial  variations  are  of  interest.  Further,  we  will  first  obtain  some  relations  involving 
the  Fourier  transforms  of  these  quantities.  Fourier  transform  (over  bar  denoting  the 

transform)  of  the  function  Fa  and  the  inverse  relation  are  given  by 

Fa  (**.  Ky>  z)= J1 C  F«  fe  y> z)  exP  G  ** x+  j  Kx  y)dx<Jy  *  (B~7) 


Fa  (x,y,z)=7TEl  Fa  (fcx,Ky,z)exp  (-j*,  x-  jar,  (B'8) 

4/r 

Similar  relations  hold  for  other  functions  and  their  transforms.  Substitution  of  equation 
(B-2)  in  (B-1 )  and  use  of  (B-7),  (B-8)  result  in  two  ordinary  differential  equations. 


-k2x-kI) Fa  =  °;  «  =  x,y  (B‘9) 

dz 

whose  solutions  are 

Fa  =  Aa  cos  k0 A  z  -  Ba  sin  kc  A  z;  «  =  x,y  (B-10) 

where  Aa  and  Ba(«  =  x>y)  are  unknown  complex  valued  functions  of  kx  and  Ky,  and  z  is 

measured  from  the  center  of  the  layer  under  consideration.  A  is  the  characteristic  root 

given  by 


such  that 


(B-1 1 ) 


Re(/l)>0  and  Tm(A)<0 

k'x,  K'r  are  nondimensional  transform  parameters,  i.e., 
k'0  =  kJ k0;  a  =  \,y 

Substituting  (B-10)  in  (B-3)  and  (B-4),  one  obtains 


(B-1 2) 

(B-1 3) 
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(B-14) 


Ex  =  "  k0  ^  (Ay  sin  k0  ^  z  +  By  cos  k0  ^  z) 
Ey  =  ko^  (AxSin  Mz  +  BxCOSko^z) 

h  x  =  J  ®  Hx  =  [(k?  -  Kx)  Fx  -  /Cx  Ky  Fy] ^ 

r* 

hy  =  j  ®  /'o  Hy  =  [(kf  -  /Cy)Fy  “  *Tx  Ky  Fx]^ 

r* 


(B-1 5) 

(B-1 6) 

(B-1 7) 


Substituting  z=+h/2  and  -h/2  in  (B-14)  one  obtains  two  expressions  which  can  be  inverted  to 

obtain  Ay  and  By  in  terms  of  eI  and  Exdefined  as 

E*  =  Ex  atz  =  ±h/2  (B-18) 

By  using  a  similar  procedure  with  equation  (B-1 5)  with  z  =  ±h/2,  Ax  and  Bx  can  be  expressed 

in  terms  of  e£  ,  where 


Ey  =  Ey  atz  =  ±h/2 

Substitution  of  these  relations  in  (B-1 0)  and  in  (B-1 6)  and  (B-1 7)  at  z  =  ±  h/2  yield 


k0A// 

Mo 


r/r-1 

Ie:1 

r  e:  ' 

\ - 

1 _ 

=  -[B] 

(kf-jcf) 

i 

1 

J*ih 

l _ 

+  KxKy 

JC 

(B-1 9) 


(B-20) 


where 


k0  *M 

Mo 


f  e:  1 

r  e:  1 

i _ 

=  [B] 

(kf-KTy) 

+  w 

rN 

1  i 

_ i 

+  /Cxxry 

[ 

1 

hx  =  hx  at  z  =  ±  h/2 


(B-21 ) 


(B-22) 


hy=hv  atz  =  ±h/2  (B-23) 

and  the  matrix  [B]  is  given  by 

cotko'ih  cosecko'ih 

(B-24) 

cosecko^h  cotko^h 

Multiplying  both  sides  of  (B-20)  by  Ky  and  those  of  (B-21)  by  -/cx ,  adding  of  the  resulting 


expressions  and  inverting,  one  obtains 
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KTyEy+KxEx  KTyhx"K*hy 

ko  _  _+  =[C] 

^ifix  ^yhx"^*hy 


(B-25) 


where 


V2uPi  ,  M  -  cosec ko^h 

[C]  =  -i^±  [B]  = J— - 

ki//0  iWo^L'00860^^^  cotko^h. 

A  similar  procedure  involving  multiplication  of  both  sides 


where 


^"^yEx 

JCxhx  +  ^yhy 

_  =[D] 

+  K*yEx_ 

JTxhx  +  *yhy 

-m'  =  JLr 

cotko^h 

X  HQX 

-  co  sec  ko'lh 

(B-26) 

s  of  (B-20)  and  (B-21)  by  k*.  and 
inversion  yields 

(B-27) 


(B-28) 
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COMPLETE  SOLUTION  FOR  THE  LAYERED  MEDIUM 

It  may  be  noted  that  relations  (B-25)  and  (B-27)  are  similar  in  form  for  all  layers  of  finite 
thickness  with  the  elements  of  the  matrices  [C]  and  [D]  differing  from  layer  to  layer  because 
of  thickness  and  material  property  differences.  At  the  interface  M+l)  of  layers  «  and  (HT), 
we  denote  the  quantities  with  a  superscript*,  i.e., 

(£)  (£+1)  (£  +  \) 

e;  =  e-  =  e: 

(B-291 

(£)  (£  +  1)  (<  +  l) 

- —  * *  * 

ha  —  ha  —  ha 

11?  ^1?  denote  the  values  of  the  quantities  at  the  top  of  layer  1  (global  z  coordinate  z  =  0 

Ea  ha 

(n  + 1)  (n  +  1) 

in  the  plane  of  the  aperture  A  in  Figure  B-1)  and  — •  and  -  are  the  values  at  the 

Ea  ha 

bottom  of  layer  n  (z =h,  +....+  h.)»  sum  of  the  thicknesses  of  layers  1  to  n. 

Now  we  add  second  of  equation  (B-25)  for  layer « to  first  of  the  equation  (B-25)  for  layer 
({+1)  for  ?  varying  from  1  to  n  and  assemble  the  following  system  of  equations  for  evaluation 

of  3  5  (/  =  1, 2, ...  w  +  l) 

Ky  hx  -  Kx  hy 


B-6 


B-7 


where 


^  (ij  =  i  ,2)  are  the  elements  of  the  matrix  [C]  (equation  (B-26)  of  layer  i).  Note 
cv 

that  the  first  of  equation  (B-25)  for  layer  1  and  second  of  (B-25)  for  layer  n  are  left  unaltered 
in  the  system  (B-3)  as  the  first  and  the  last  equation. 

A  similar  procedure  (addition  of  second  of  equation  (B-27)  for  layer  l  to  first  of  (B-27) 


for  layer  (C+l)  yields  the  following  system 


of  equations  for  evaluation  of. 


<0  <0 

Kx  hx  +  Ky  hx 


(i=  1,2,...,  n+1). 


(B  -  32) 
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0) 

*12 

(1) 

1^*22 +*11, 

(2) 

d\2 


(2) 


12 

'© 

Kd22  +dlV 


0 


(3) 

f12 


(/-l)  f(/-l)  (/)' 


*12 


l  *22  +  dlV 


0 


*12 

'(»-0  (»r 

0-0 

w 

*12 

<  *12  +  *11, 
1/21 

*12 

(») 

VV 

*12  • 

YV 

*22. 

(B-33) 


Note  that  ^  *  (ij-1,2)  are  the  elements  of  the  matrix  [D]  (equation  (B-28)  for  layer  i).  Again, 
da 

first  of  equation  (B-27)  for  layer  1  and  second  of  (B-27)  for  layer  n  are  left  unaltered  as  the 
first  and  last  equation  in  the  system  (B-32). 

CASE  1.  LAYER  n  BACKED  BY  A  CONDUCTOR 


For  a  layered  system  with  the  n,h  layer  backed  by  a  conductor  (Figure  B-1),  Ex  and  Ey  are 
identically  zero  at  the  back  of  layer  n  (global  z  =  h,  + ...  hn)  and,  therefore,  the  last  terms  in 
the  vectors  on  the  left  hand  sides  (B-30)  and  (B-32)  are  zero.  Therefore,  each  term  of  the 
vector  involving  the  transforms  of  the  magnetic  fields  tangential  to  the  interfaces  (or  rather 
their  derivatives)  appearing  on  the  right  hand  sides  of  (B-30)  and  (B-32)  can  be  evaluated  in 
terms  of  the  transforms  of  tangential  electric  fields  at  the  plane  of  the  ap&ure. 

In  particular,  the  transforms  of  the  tangential  magnetic  fields  at  the  plane  of  the  aperture 


are 


(i)  0) 

— *  — * 

Kv  hx  -  Kx  by 


0) 

A,  KyEy  +KxEXy 


(B-34) 


and 


B-9 


0)  0)  £*(  0)  0) 

ATX  hx  —  Ky  hy  =  A2  j Ey—KyEx 


where 


Ai  =  Det 


r(G)i  r(Gj 

5i  =  Adjoint  of  in  determinant  of 


A,  =  Det 


m  m 

52  =  Adjoint  of  in  determinant  of 


(B-35) 


(B-36) 


CASE  2.  LAYER  n  BACKED  BY  A  HALF  SPACE 


If  the  layered  system  is  not  backed  by  a  conductor,  but  by  a  half  space  of  material  me¬ 
dium  (n+1),  then  the  tangential  electric  fields  are  not  zero  at  the  back  of  layer  n.  However, 
the  transforms  of  electric  fields  are  related  to  those  of  the  magnetic  fields  at  the  interface 
(n+1)  between  layer  n  and  medium  (n+1)  by  the  following  relations. 


f  ("i1)  ("+0  ("_+ 0 

T_KV  Ey  ~Kx  E*  J  IS  hx  -kx  hy 


(B-37) 


(w  +  l)  (w  +  l)"l_  (w  +  l)  (n  +  l) 


°l  —  ATV  Ey  +  fc  Ex 


K+\Eo  \Kx  hx  +Ky  hy 


(B-38) 


Subscripts  (n+1)  indicate  quantities  for  the  medium  n+1.  Substituting  these  relations  in  the 
left  hand  sides  of  the  last  equations  in  the  systems  (B-30)  and  (B-32)  and  transferring  them 
to  the  right  hand  sides  of  these  equations,  the  last  terms  in  the  vectors  on  the  left  hand 


sides  become  zero.  However,  the  terms 


are  changed  to 


'(n+lX«+!) 


/(n+lX«4l) 


"(G) 

of  the  matrices  and 

C 
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(B-39) 


(G) 

C(n+1Xn+l)  Crj  ^n+1  Mo 


(G) 


C22 

(n) 

=  +• 


JiVi 


(B-40) 


d(n+lXn+l)  d22  ■^n+1  ^0 

After  this  modification,  inversion  of  (B-30)  and  (B-32)  can  be  done  in  the  same  manner  as 
discussed  for  case  1  (see  equation  (B-34),  (B-35)  and  (B-36),  which  apply  to  this  case  too, 

with  the  modifications  of  j  and  |  ). 


r<°>!  and 

[(G)] 

u 
— 1 

D 
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TRANSFORM  INVERSIONS  AND  SOLUTION  PROCEDURE 


From  equations  (B-34)  and  (B-35),  it  follows  that 


(B-41 ) 


(B-42) 


These  transforms  can  be  inverted  using  the  relation  (8)  to  evaluate  ^  and  ,  provided  ^ 

fix  hy  Ex 


and  ,  are  known  over  the  aperture  (they  are  zero  elsewhere).  For  an  approximate  solu- 

Ey 


tion,  the  tangential  electric  fields  at  the  aperture  may  be  assumed  to  be  the  sum  of  the  inci¬ 
dent  and  reflected  parts  of  the  dominant  TE  mode  for  the  wave  guide  or  the  co-axial  cable. 
The  corresponding  magnetic  field  can  then  be  evaluated  for  calculation  of  the  aperture  ad¬ 
mittance,  which  can  be  used  for  suggesting  guidelines  for  nondestructive  examination.  In 
the  following  sections,  we  will  discuss  procedures  for  such  calculations.  It  should  be  noted, 
however,  that  for  obtaining  a  correct  solution  other  modes  in  the  wave  guide  have  to  be 
considered  and  all  the  tangential  fields  at  the  aperture  should  be  matched  with  those  of  the 
wave  guide.  For  this  purpose,  it  is  necessary  to  express  the  tangential  fields  at  the  aperture 
in  the  wave  guide  atz  =  0  as 

E (x,  y, o)=  Vio  e0  (x,  y)+  £  Vm  en (x,  y)  (B-43) 

n=0 


H  (x,  y,  o)=  Y0  Vio  h„  (*,  y)  ■ -  £  Yn  V™  hn  (x,  y) 

n=0 


(B-44) 


Vio ,  Vro  are  the  magnitudes  of  the  dominant  incident  and  reflected  modes,  Vm  are  the  mag¬ 
nitudes  of  other  reflected  modes.  Y„  is  the  admittance  of  the  wave  guide  in  the  11th  mode, 
which  is  known.  e„  and  hn  are  the  orthonormal  mode  vectors.  The  normalized  aperture 
admittance  for  the  dominant  mode  is  given 
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(B-45) 


Y  _  Vjo  ~  Vro 
Yo  Vio  +  Vro 


For  numerical  solution,  a  finite  number  of  terms  (n  <  N)  are  retained  in  (B-43)  and  (B-44)  and 
Vm  (n=0,l...N)  are  expressed  in  terms  of  Vio  by  matching  the  magnetic  fields  obtained  from 
(41)  and  (B-42)  with  those  from  (B-44)  by  collocation  (point  matching  over  a  chosen  number 
of  points)  or  other  error  minimization  procedures.  It  may  be  noted  that  some  of  the  fields 
may  be  singular  at  the  aperture  boundaries  when  the  loss  tangents  for  the  material  dielectric 
constants  are  zero  (see  Reference  2)  and  other  improved  techniques  (singular  integral  equa¬ 
tions)  may  be  required  for  solution.  When  the  loss  tangents  are  not  small,  one  term  solution 
(n=0)  may  be  adequate  (see  Reference  1)  and  the  procedure  often  yields  good  correlation 


with  experimental  data.  For  such  solution,  one  may  choose 
Vio  +  V^l 

vm=0  ;  «>0 


(B-46) 

(B-47) 


and  it  is  only  necessary  to  find  Vio  -  Vro  from  the  dominant  magnetic  field  (see  eqn  (B-44)) 
and  use  (45)  to  find  the  normalized  aperture  admittance.  The  case  of  a  rectangular  wave 
guide  and  that  of  a  co-axial  cable  are  considered  in  the  next  two  sections. 
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APPROXIMATE  SOLUTION  FOR  RECTANGULAR  WAVE  GUIDE 


The  wave  guide  is  assumed  to  be  of  lengths  a  and  b  in  x  and  y  directions,  respectively, 
(Figure  B-2)  with  a>b.  Considering  only  the  dominant  mode  in  the  guide,  which  is  the  TE10 
mode,  and  neglecting  all  other  reflected  modes,  one  may  choose 


0) 

^  JVC 

Ev  —  Cos — 
y  a 


(1) 

Z= o 


(B-48) 


-  a/2  <  x  <  a/2  and  -  b/2<  y  <  b/2 


(B-49) 


Figure  B-2.  Rectangular  Wave  Guide 


Note  that 

^  ^  (B-50) 

Ex  =  Ey  -0,\x\>  a/2  or  [y|  >  b/ 2 

because  of  the  presence  of  the  conducting  flange  (Figure  B-1).  Using  equation  (B-7)  and 

(1) 

(B-48  -  B-50)  as  well  as  the  condition  that  *  is  an  even  function  of  x  and  y 

jb„ 
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0) 

%=\-a2n  (In  cosnx/acosKxxcosKyydxdy 
An  co$Kxal2sinKybl2 


(B-51) 


a  (;r2/a2-jc2)t. 


0) 

£*x=0 


<B-52) 


0) 

//* 


Use  of  equations  (B-52),  (B-41),  and  (B-8)  and  the  condition  that  -/  is  an  even  function  of 
kx  and  Ky  yields 

(1) 

v  J  (B-53) 


(1)  (1)  ka  ^  „  1  U  2  S2  2 

jw. EL  C.  (v>+^)[4|  Kr +  Aj  *' 


Ey  cosxrxx  cosxrJj'  dKxdtc 


Now  we  note  that  within  the  aperture  can  be  expanded  in  the  following  form  (as  dis- 
H 


cussed  in  the  previous  section) 


(1)  0)  (1) 

tf*=/C«>s*r/<r+2  H*xJ„(x,y) 


(B-54) 


where  fn(x,y)  are  the  orthogonal  mode  functions  for  higher  modes  of  the  wave  guide. 
H*c an  be  evaluated  as 


0)  9  0) 

H'x°=—\°ai2  t'bu  H*xcosnx/adxdy 


(B-55) 


Substitution  of  equation  (B-51)  in  (B-53)  and  the  resulting  expression  in  (B-55)  as  well  as 
evaluation  of  the  integrals  with  respect  tox  and  y  (see  (B-55))  yields 


(1)  8*  a 

Z.IdK,dK, 


(B-56a) 


where 
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(B-56b) 


1  = 


1 

~Si 

Wx+jcj) 

_Ai 

I  2  4.  2 

'Ky  + Kx 

A2  . 


cos2  Kx  a/  2  sin2  xry  b/  2 


^-xWj 


2 

Ky 


Noting  that  the  aperture  admittance  Y  is  the  ratio  of  H*m  and  (the  latter  is  unity  as  de¬ 
fined  by  equation  (B-48),  see  equations  (B-43-45)  and  accompanying  discussions),  one  ob¬ 


(B-57) 


tains 

(1) 

Y  =  G+jB=H'xo 

which  is  given  by  (B-56a,  b).  Y0<  the  admittance  of  the  wave  guide  for  the  dominant  mode 
is  given  by 

Yo  =  -A/k2-^2/a  2/®/*0  (B'58) 

ko,  Po  are  the  wave  number  (propagation  constant)  and  the  permittivity  of  the  free  space. 
The  normalized  admittance  of  the  aperture  is,  therefore,  given  by 

0) 

Y  H*xo0)juo 


y  =  g  =  ib  = 


yjlc*  -7T2  la7 


(B-59) 


0)  . 


H 


\  is  given  by  (B-56).  Other  quantities  of  interest  are  given  below. 


Reflection  Coefficient 

Rc  =  (l-y)/(l+y)=€'ei*  (B-60) 

<|)  is  the  phase  angle 
T  is  the  magnitude 

Voltage  Standing  Wave  Ratio  (dB) 

VSWR  =  20  logo  (  (  1  -O  /  (  1  +1")  )  (B-61) 

The  integral  in  (B-56a)  can  be  converted  to  a  sum  of  two  infinite  series  by  the  use  of  calcu¬ 
lus  of  residues.  The  residues  have  to  be  evaluated  at  the  zeroes  of  Ai  and  A2.  This  proce¬ 
dure  is  imperative  when  the  materials  have  zero  loss  tangents.  If  the  loss  tangents  are  not 
zero,  the  integral  can  be  evaluated  numerically.  For  this  purpose,  the  following  transforma¬ 
tions  may  be  utilized  along  with  the  property  that  the  integrand  is  even  in  kx  and  Ky. 
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jCx  =  ko^cos(9 


0<X-<oo 


Q<6<n!2 

«ry  =  ko^sin^ 

k  =  i 7R/(1-R)  ;0<R^1 


(B-62) 


77 = a  chosen  multiplier 

The  integral  can  then  be  written  as  follows.  The  last  expression  uses  a  quadrature  formula, 
w,  and  w'M  being  the  weights. 

r»r.l(>Cx,rry)d^d^y  =  4k^j:/2l(/rx^yVd^d/r 

=  4koJoJo/2/(ko^cos0,ko^sin  (B-63) 


=  4ko£w/  Z  w'  l(k0K-/ cos0m,koKt sin Om),**  ^ 
/=J  m=J  M'Rr/ 

_  *7R /  0  - 1 


where  tct  ~  ~  ^-lv  -L 
1-R/ 

If  one  uses  trapezoidal  rule  with  equidistant  integration  points  in  the  domains  (0  4Rs  1 
and  0  <  0  <  nil),  then 


r,  =  /?/(L+1)  ;  ^  =  1,2,...L 

_(m-l);r  ;  m  =  l,2,....M 

6>m_2(M-l) 

W/=(L+i) 

w'f  =  ;r/4(M-l)  m  =  1,M 
=  tt/2(M-1)  2  <  m  <  M- 1 


(B-64) 


It  may  be  noted  that  the  points  R  =  0  and  1  (k  —  0  and  °°)  are  omitted  in  the  quadrature  for¬ 
mula,  since  the  value  of  the  integrand  is  zero  at  these  points. 
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The  expression  I  in  (B-63)  can  be  written  in  the  following  form  in  terms  nondimensional- 


ized  quantities 


,  2t _  S\  ; 
kol - 7~C0S 

A, 


1*.  62  .  cos2 (* a'cos 6 / 2)sin2 (k b'sin 0/2) 

\  0  + — sm  0  - 7— — ”  1 


(x2  -  k  a'2  cos2  oJk2 sin2  0 


(B-65) 


a'  =  k0a,  b'  —  kQb  (B-66) 

Closed  form  expressions  of  8a  and  A*  (a  =  1,2)  appearing  in  56(b)  and  (65)  are  given  below 
for  a  few  cases.  The  following  notations  are  used  in  the  expressions  (subscript  denotes  the 
layer  or  backing  medium  p). 


p  =  — — -  for  a  =  1 


(B-67) 


for  a  =  2 


(B-68) 


Cp  =  cos(k0Aphp)  (B-69) 

Sp  -  sin(k0AJ^rJ>)  B‘70) 

The  expressions  which  follow  are  applicable  to  a  =  1  and  2,  but  appropriate  values  of  pp 
(B-67)  for  a  =  1  and  (B-68)  fora  =  2  have  to  be  used. 


A.  Single  Layer  (p=l)  Backed  by  a  Half  Space  (p=2) 

Sa  =  Px  G+j&S.  (B7 

Aa  =  -^S,  +  jAAC. 

For  a  single  layer  with  conductor  backing,  the  terms  with  j  p2  have  to  be  omitted. 


(B-71a,b) 


B.  Two  Layers  (p=l,2)  Backed  by  a  Half  Space  (p=3) 

Sa  =  px  p2a  C2 -/?2s,s2+j  &(0iC.  S2  +  £2S,  C2) 

Aa  =  -^,^2tos,C2  +  /?2C,S2) 

+j/M.(-A>s,s2+£2cIc2) 

For  two  layers  with  conductor  backing,  the  terms  with  jp3  have  to  be  omitted. 


(B-72a,b) 


B-18 


C.  Three  Layers  (p=l,2,3)  Backed  by  a  Half  Space  (p-4) 

Sa  =  P3(pxP 2 C, C2 C3 -  A P3 C1S2S3- Pi S, S2 C3 
-P2P3SiC2S^  +  jP4(PlP2CiC2S3 
+PlP3CiS2C3-P22SlS2S3+P2P3SxC2C3) 

(B-73a,b) 

Aa  =  P\  P3  {P\  P3  Si  S2  S3  ”  P\  P2  Si  C2  C3 -  P\  Ci  S2  C3 
- P2P3O  C2S3)+  j P4PX  (-  Px  p2 S.  C2 S3 
-  Px  p3  S.  S2  C3  -  Pi  C,  S2  S3 +  P  2  P 3  Cl  C2  C3) 

For  three  layers  with  conductor  backing,  the  terms  with  jp4  have  to  be  omitted. 

D.  Four  Layers  (p-1,2,3,4)  Backed  by  a  Half  Space  (p-4) 

Sa = P4\P2  teCA — P2P ~  Pa^ 3^ a) 

+ a^aXMc. + AC,s.) ) 

+  JPslfil  iP &P2  —  Pl^Pl X/^3^3^4  +  PP3P4) 

+  AfoCA  +a%Xac#.-AAA)  ] 

(B-74a,b) 

Aa  =  Pi  iPP\^i  +  Pi^Pi HP3C3C4  ~~ P4S3S4) 

+ a(AAA  -  fi2c,cMs£, + aca)  ] 

+JAAI-A  (Me,  +AAC.XACA + AAQ) 

-  -  PAcMcfi,  -  0,SA)  ] 

For  four  layers  with  conductor  backing  the  terms  with  jj34  have  to  be  omitted. 

It  may  be  noted  that  the  expressions  for  5  and  A  for  more  than  four  layers  (without  or 
with  a  half  space  backing)  are  not  incorporated  in  the  code  written  for  computation  pur¬ 
poses  and,  hence,  the  code  cannot  handle  greater  than  four  layers.  Similar  expressions  can 
be  easily  obtained  and  inserted  in  the  code,  if  desired. 
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SERIES  SOLUTION  FOR  A  RECTANGULAR  LAYER  BACKED  BY  A  CONDUCTIVE  SHEET 


As  mentioned  earlier,  the  integral  (B-53)  can  be  expressed  as  a  sum  of  two  series  by  the 
use  of  the  calculus  of  residues.  The  residues  have  to  be  calculated  at  the  complex  zeroes  of 
A)  and  A2.  For  the  case  of  a  single  layer  of  thickness  hi,  the  zeroes  of  Ai  are  of  interest, 
which  are  the  roots  of 

sin(MA)=°  (B'75) 

and  the  series  solution  for  normalized  guide  admittance  is  given  by 

jX 

■Jl-n  /a'2  L  °  '  bh'  n~0-' 


y ,  = 


.....  4^1  -  a„  exp(-  jk'nb '  /  2)sin  (jk'„b'/  2) 

-cot k0h,  +  —  Z  - £75 


(B-76) 


where 


sr 

j 

k.i'.A’l 

k, 

On 


=  \sr-n2  la'1  —n2n2  Ihf]'2  with  negative  imaginary  part;  n  =  0,1  ,... 

=  complex  dielectric  constant  of  the  lossy  plate 

=  yH  (B-77) 

=  [a,b,hi]  ko,  a,b  are  the  wave  guide  dimensions  and  hi  is  the  plate  thickness 

=  free  space  wave  number 
=  1/2  for  n=0  and  1  forn>0 


The  series  and  the  integral  have  been  found  to  give  identical  results  and  are  compared  in 
Appendix  C  with  some  results  available  in  literature  for  the  single  layer  case. 
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SOLUTION  FOR  A  CO-AXIAL  CABLE 


For  a  co-axial  cable,  the  aperture  is  in  the  form  of  an  annular  region  bounded  by  two 
concentric  circles  of  inner  and  outer  radii  a  and  b,  respectively  (Figure  B-3).  In  this  case,  the 
geometry  is  referred  to  a  cylindrical  coordinate  system  p,y  (See  Figure  B-3)  and  z,  where 

x  =  pcosy/  (B-78) 

y  =  psiny/ 

and  we  restrict  our  attention  to  the  case  of  an  axisymmetric  electric  field  incident  on  the 
aperture.  Thus,  in  the  layered  medium 

Ex~EPip,  z)cosy/  .  (B-79) 

Ey  =  Ep(p,z)sMV' 

and 

Hx  =  -n,{p,^¥  (B-80) 

Hy  =  Ky,{p,2)<x>sV 


Figure  B-3.  Co-Axial  Cable 

Using  (B-78)  and  the  first  two  equations  of  (B-62)  expressing  the  transform  variables  kx  and 
Ky  in  terms  of  (k,0),  one  can  write 
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exp(j  ACx  x  +  j  Ky  y)=  exp(j  kn  tc o  cos(^  -  0)) 


(B-81 ) 


The  right  hand  side  of  (B-81)  can  be  expressed  in  a  series  form  given  below  [3] 

expO  ko  *P  cos<y  -  0))  =  Jo  (ko  *p)+  2  Z  0  )”  Jn  (k0  xp)cosn(yr  -  0 )  (B-82) 

n=l 

Substituting  (B-82)  in  (B-7)  and  changing  the  integration,  with  respect  to  x  and  y  to  that 
with  respect  to  p,  y,  one  can  evaluate  the  Fourier  transforms  of  the  quantities  in  (B-79)  and 
(B-80).  Performing  the  integration  with  respect  toy,  the  following  results  are  obtained. 

^  =  2n)Epcos0 
F„  =  27rj  E„  sin# 

_y  *  (B-83) 

=  -  2;r  j  sin  0 

Hy  =  2;r  j  Hv  cos  <9 

where  ~  over  a  quantity  indicates  its  Hankel  transform,  i.e., 

E/,  =  JT^EPJi(koK'p)d/?  (B-84) 

Ji  denoting  the  Bessel  function  of  the  first  kind. 

It  follows  form  (B-83)  that  because  of  axisymmetry 
xrxEy-KyEx  =  0 

and  (B-85) 

Therefore,  one  needs  to  use  only  (B-30)  to  solve  the  problem.  Equation  (B-32)  is  identically 
satisfied.  Using  (B-34)  for  the  case  with  conductor  backing,  or  its  modified  form  for  the  half 
space  backed  case  (as  discussed  in  case  2  earlier),  one  obtains  the  following  relation  be¬ 
tween  the  Hankel  transforms  of  the  magnetic  and  electric  fields  on  the  top  surface  of  the 
layered  medium. 


(1)  S'  5> 

K=jWoK  =~~txk°E*p 


(B-86) 


As  discussed  earlier,  a  correct  solution  can  be  obtained  by  assuming  the  electric  field  at  the 
aperture  as  a  sum  of  the  incident  and  reflected  parts  of  all  possible  modes  in  the  co-axial 
cable  (the  field  is  zero  outside  the  aperture),  and  matching  the  magnetic  fields  at  the  aper- 
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ture  with  those  obtained  by  inverting  (B-82)  (see  12]).  For  obtaining  an  approximate  solu¬ 
tion,  when  the  loss  tangents  of  the  layers  are  not  zero  (or  not  too  small),  one  can  consider 
only  the  dominant  incident  and  reflected  modes  in  the  co-axial  cable,  i.e.. 


0) 

K=l/p 

=  0 


a<  p<b 
p<a,p>b 


(B-87) 


Substituting  (B-87)  in  (B-84)  and  performing  the  integration  between  p=a  and  p=b,  one  ob¬ 


tains 


0) 

K = -p.m: kwv 


(B-88) 


(l) 

Use  of  this  result  in  (B-86)  and  inversion  of  the  Hankel  transform  in  (B-86)  result  in 

Hv 


(1)  , 

Hl  =  \:  XKKb)  -  J0{K^)V^P)d{K^) 

“l 


(B-89) 


JG>Po“v 

Now  we  note  that  the  magnetic  field  in  the  wave  guide  at  the  plane  of  the  aperture  can  be 
expanded  in  the  form 


(i)  0)  „  (i) 

K=h'Jp+^  Kw 

«=1 


(B-90) 


where  fn  (p)  are  the  orthogonal  mode  functions  of  higher  order,  satisfying  the  condition 

£f„(p)<ip=o  (B-91) 

Substituting  (B-89)  in  (B-90),  integrating  with  respect  to  p  (from  a  to  b)  and  using  (B-91), 
one  obtains 


0)  1  *  0) 
H*°-  C,n{b I a)^a  H*dp 

= — - — r 7  a 

jcop/n(b/a) 

where 


(B-92) 
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(B-93) 


&\K 

a'  =  k0aandb'  =  k0b 

As  discussed  in  the  case  of  the  rectangular  wave  guide  problem,  the  aperture  admittance, 

Y,  is  the  ratio  of  ^  and  £!?  (the  latter  is  chosen  as  unity  in  (B-87))  and,  therefore, 

E„„ 


(1) 

Y=K 


(B-94) 


yo 


which  is  given  by  (B-92,  B-93).  Y0,  the  admittance  of  the  co-axial  cable  in  the  dominant 
mode,  is  given  by 


Y„  =  G  +  jB  =  — =koJv£ 

0)flc  0)flc 


(B-95) 


where  the  quantities  with  subscript  c  indicate  properties  of  the  core  of  the  co-axial  cable. 
Therefore,  the  normalized  admittance  can  be  written  as 


y  ~  g + jb = ~~  =  ^I,dy 

Yc  Mofjl\h/ahl£c 


(B-96) 


Other  quantities  of  interest,  such  as  reflection  coefficient  R*  (y  =  (1-R<)  /  (1+Rc)),  etc.,  can  be 
computed  in  the  same  manner  as  discussed  for  the  rectangular  wave  guide  problem. 
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COMPUTER  CODE 


A  computer  code  ADMIT  has  been  written  to  evaluate  the  integrals  and  results  for  the 
case  of  one  term  of  approximation  and  non-zero  loss  tangents.  The  following  input  proper¬ 
ties  are  fed  in  NAMELIST/INPUT. 


NF 

FTNT 

DF 

EPSR  (1),  EPSI  (I) 

H(I) 

A,B 

ITYP 

NMAT 

NLAY 


CSPR,  CSPI 

NZET 


ETZET 

MPHI 


=  No.  of  Frequencies 
=  Starting  Value  of  Frequency  in  GHZ 
=  Increment  of  Frequency  in  GHZ 

=  Real  and  Imaginary  Part  (Negative)  of  Dielectric  Constant 
for  Material  1=1,  NMAT 
=  Thickness  of  Material  Layer  I 

=  Dimension  of  Rectangular  Wave  Guide  (A  >  B)  or  Inner  and 
Outer  Radii  of  Co-Axial  Cable  (A<B)  in  Inches 
=  1  for  Rectangular  Wave  Guide,  =  2  for  Co-Axial  Cable 
=  Total  No.  of  Materials 

=  No.  of  Layers  of  Finite  Thickness  (limited  to  <4  at  this 
time) 

(Note  that  NLAY =NAMT  for  Conductor  Backed  Case  and 
NALAY = NMAT-1  for  Half-Space  Backed  Problem) 

=  Real  and  Imaginary  Parts  of  Dielectric  Constants  of  Core 
for  Co-Axial  Cable  Case 

=  No.  of  Integration  Points  for  Integration  (B-63)  with 
Respect  to  R,  Value  of  L  in  last  of  equation  (B-63).  Note 
that  a  similar  transformation  and  quadrative  formula  is 
used  for  evaluation  of  the  integral  in  (B-92)  for  co-axial 
cable. 

=  The  multiplier  used  in  (B-62) 

=  No.  of  Integration  Points  for  integration  (B-63)  with 
respect  6,  value  of  M  in  second  of  equation  (B-64) 


ll  * 


NMH 

NHE 

DH 

DEPSR,  DEPSI 
UN,  IDDT 

KOPT 

rxx(I),ryy(I);I=i,nf 


RXZ(I),RYZ(I),I=1,NF 


=  Sequence  No.  of  Material  Layer  whose  thickness  or 
dielectric  constants  are  to  be  varied,  if  desired 

=  No.  of  Thickness  Increments  of  Material  Layer  NMH  as 
desired 

=  Thickness  Increment 

Note  that  frequency  is  varied  (N=l, ...  NF)  first  for  each  set 
of  layer  thicknesses  or  properties  (NH=1,  ...NHE). 

=  Increments  in  Dielectric  Constants  of  Material  NMH  as 
desired 

=  Print  Options  for  checking  purposes.  No  input  is  required 
if  such  printings  are  not  needed.  They  are  set  to  zero. 

=  For  calculation  of  changes  in  the  lift-off  and  orthogonal 
directions  (see  Appendix  E) 

=  Real  and  Imaginary  Parts  of  Reflection  Coefficients  for  No 
Defect 

=  Real  and  Imaginary  Parts  of  Reflection  Coefficients  for  a 
Very  Small  Defect  Thickness 


All  inputs  are  straightforward  except  those  required  for  integration.  A  value  of 
ETZET  =  4.0  should  be  adequate  for  all  cases.  It  has  been  found  that  MEES  of  the  order  30 
or  40  is  required  for  convergence  of  the  integral  with  respect  to  0.  Note  that  although  this 
integration  (with  respect  to  0)  does  not  appear  in  the  case  of  co-axial  cable  (eqn  (88)),  the 
Bessel  functions  are  evaluated  in  the  code  using  an  integration  scheme.  MPHI  =  30  to  40 
was  found  to  be  all  right  for  such  evaluation.  The  value  of  NZET  is  very  crucial  and  it 
should  be  of  the  order  of  12000  when  the  loss  tangents  are  small  (of  the  order  of  -0.001). 
For  loss  tangents  of  the  order  of  -0.1,  NZETA  -  1000  may  be  adequate.  NZETA  =  6000 
may  work  for  a  loss  tangent  of  -0.01.  It  should  be  noted  that  four  integrals  are  evaluated 
for  integration  points  (NZET,  MPHI),  (NZET,  MPHI/2),  (NZET/2,  MPHI)  and  (NZET/2,  MPHI/2) 
and  printed  for  checking  convergence. 

Output  consists  of  input  data,  normalized  and  actual  admittances,  complex  reflection 
coefficient,  its  magnitude  and  phase  (in  degrees)  and  VSWR.  Plot  files  are  created  to  give 
real  and  imaginary  points  of  admittance  (normalized  or  actual)  as  well  as  magnitude  and 
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phase  of  reflection  coefficient  versus  frequencies  or  thickness  variations  for  a  particular 
layer.  The  plot  file  ADMAT.PLTA  contains  frequencies  (or  thicknesses  of  material  NMH  if 
NF  =  1),  real  and  imaginary  parts  of  normalized  admittances  and  real  and  imaginary  parts  of 
actual  admittances  in  five  successive  columns.  The  plot  file  ADMIT.PLTB  has  frequencies 
(or  thicknesses  of  material  NMH  if  NF  =  1),  real  and  imaginary  parts  of  the  reflection  coeffi¬ 
cients  and  VSWR  (db)  and  phase  (degrees)  in  five  columns.  The  plot  file  ADMIT.PLTC  con¬ 
tains  frequencies,  the  projections  of  the  reflection  coefficient  vectors  for  no  defect  case  on 
to  the  lift  off  vector  and  the  orthogonal  component  as  well  as  the  changes  in  the  projections 
when  there  is  a  defect  (from  no  defect  case).  The  file  ADMIT.PLTD  has  frequencies  and  the 

direction  cosines  of  the  lift  off  vector  and  those  of  the  orthogonal  component  in  five  col- 

# 

umns. 

A  list  of  the  code  follows. 
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LISTING  OF  THE  CODE  ADMIT 


C  PROGRAM  ADMIT 

C  COMPUTES  ADMITTANCES  OF  RECTANGULAR  WAVE  GUIDE  OR 
C  CO-AXIAL  CABLES  WITH  FLANGES  FACING  A  LAYERED  MEDIUM 


C 

C  PC  COMPILE  INCREASE  STACK  SIZE  TO  HEX  COO  =  3K  decimal 
C  ACTUAL  COMPILE  LINE  BELOW:  MS-FORTRAN  VERSION  5.1 
CFLFPi/F  COO  ADMIT.FOR 

C  PC  CHANGE ADMIT.PLT A,  ADMTT.PLTB  TO  ADMIT.PTA,  ADMIT.PTB 
C  DUE  TO  FILE  EXTENSION  LIMITATION  OF  3  CHARACTERS 
C 


Q********************************************************************* 


C  INPUT  DATA 


c ******************************************************************** 


C  ALL  DIMENSIONSA,B  AND  H  IN  INCHES,  FREQUENCIES  IN3HZ 

C  A,B=RECT ANGULAR  WAVE  GUIDE  DIMENSIONS  A  LARGER 

C  A,B=INNER  AND  OUTER  RADII  OF  CO-AXIAL  CABLE 

CNF=NO  OF  FRQUECIES,  FINT=  STARTING  VALUEGHZ,  DF=  INCREMENT  GHZ 

C  EPSR(I),  EPSI(I)=  REAL  AND  NEGATIVE  IMAGINARY  PARTS  OF  DIELECTRIC 

C  CONSTANTS  OF  MATERIAL  R=1,NMAT 

C  CSPR,  CSPI=  REAL  AND  IMAGINARY  PARTS  FOR  DIELECTRIC 

C  CONSTANT  FOR  CORE  OF  CO-AX 

C  ITYP=1  FORRECTAGULAR  GUIDE,  =2  FOR  CO-AX 

C  H(I)=THICKNESS  OF  MATERIAL  R=1,NMAT,  ANY  LARGE  NO  FOR  HALF  SPACE 
C  NLAY=NO  OF  LAYERS, NMAT=NO  OFMTLS 

C  NMH=MTL  NO  WHOSE  THICKNESS  OR  PROPERTIES  TO  BE  VARIED 
CNHE=NO  OF  THICKNESS  OR  PROPERTY  INCREMENTS  OF  TfflMTL  AND 


C  DH,DEPSR.DEPSI=VALUE  OF  INCREMENTS  AT  EACH  STEP  FOR  ABOVE 
C  NLAY=NMAT  FOR  CONDUCTOR  BACKED  CASE,  NMAT-1  FOR  HALF  SPACE  BACK 


C  NZET=  NO  OF  INTERVALS  FOR  INFINITEZETA  INTEGRATION 


C  ETZET=  P  ARAMETERJZETA=ETZET*R/(  1  -R),  0<R<1 

C  MPHI=  NO  OF  INTEGRATION  INTERVALS  FOR  PHI  INTEGRATION  0<PHI<PI/2 
C  BOTHNZET  ANDMPHI  MUST  BE  EVEN  INTEGERS 


C 

C  ITN=1  FORPRINTINTING  VALUE  OF  LAMBDA  AND  T AN( AMBD A*H)-CHECK 
C  IDDT=1  FOR  PRINTING  FUNCTIONS  FOR  CHECK 


Q************************************* 

C  OUTPUT  ADMIT.OUT  IS  SELF  EXPLANATORY 

q************************************************************************ 

C  PLOT  FILE  ADMIT.  PTA  CONTAINS  FREQUENCIES  (ORHICNESSES  OFMTLNMH  IF 
CNF=1)  AND  CORRESPONDING  NORMALIZEDCOMPLIANCES  AND  ACTUALCOMPLIANCES 

Q************************************************************************ 

C  PLOT  FILEPTB  CONTAINS  FREQUENCIES  (OR  THICKNESS  OIMTLNMH  IF 
CNF=1)  AND  CORRESPONDINGRFLECTIONCOEFFS,  VSWR(DB)  AND  PHASE(DEG) 
£*********************************************************************** 


IMPLICITREAL*8  (AH.O-Z) 
IMPLICfflNTEGER*4  (I-N) 
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CX)MPLEX*16AMB(10),AD(2,10)3MY,CPY,TN(10),TOP,BOT,TOP1, 

lBOTl,TOP2,BOT2,DDT(2)5UMl,SUM2,ARGC,XCl,XC2,ZSUM(4)5MULT, 

1EPSC(10),CPSC 

C 

DIMENSIONEPSR(10)£PSI(10),H(10),HNON(10),EXRA(10)EXIA(10) 

C 

NAMELIST/INPUTNF.FINT,DF,EPSR,EPSI,H,A,B,NZET,MPHI,ETZET,ITYP, 

1NMAT,NLAY,CSPR,CSPLNMH,NHE,DH,DEPSR,DEPSI,ITN,IDDT 

C 

OPEN(UNIT=5,FILE=,ADMIT.DAr,STATUS='OLD') 

OPEN(UNrr=6,FTLE='ADMIT.OUT,STATUS='UNKNOWN') 

0PEN(UNIT=7,FILE=,ADMIT.PTA',STATUS='UNKN0WN') 

C  PLOT  FILE  PTA  CONTAINS  FREQUENCIES  (OFTHICNESSES  OFMTLNMH  IF 
CNF=1)  AND  CORRESPONDING  NORMALIZEDCOMPLIANCES  AND  ACTUA1COMPLIANCES 
OPEN (UNIT =8  ,FILE=* ADMIT.PTB',  ST ATUS='UNKNO  WN') 

C  PLOT  FILEPTB  CONTAINS  FREQUENCIES  (OR  THICKNESS  OIMTLNMH  IF 
CNF=1)  AND  CORRESPONDINGRFLECTIONCOEFFS,  VSWR(DB)  AND  PHASE<DEG) 

NHE=1 

ITN=0 

IDDT=0 

DH=0.0D0 

DEPSR=0.0D0 

DEPSI=0.0D0 

RE  AD^, INPUT) 

IF  flTYP.EQ.  1)  GO  TO  1 

WRITE(6,609)CSPR,CSPI 

CALLSQROOT(CSPR,CSPI) 

CPSC=DCMPLX(pSPR,CSPI) 

1PI=3.1415926D0 

PI2=PI*PI 

SMULT=DCMPLX0.ODO, 1 .0D0) 

KMAX=2 

IF<JTYP.EQ.2)  KMAX=  1 

AN=NZET 

ANI=1.0D0/AN 

AM=MPHI 

AMI=0.5D0*PI/AM 

NZET1=NZET-1 

MPHU=MPHI+I 

WRITE(6,601)NF.FINT,DF 

WRITE(6,602)A,B 

WRITE(6,603)NZET,MPHT,ETZET 

IFdTYP.EQ.l)  WRITE(6.604)ITYP 

IF(ITYP.EQ.2)  WRITE(6,605JTYP 

WRITE(6.606)NMAT.NLAY 

IF(NHE.EQ.  1)  GO  TO  2 

WRITE(6.620)NHE 

H(NMH)=H(NMH)-DH 

EPSR(NMH)=EPSR(NMH)-DEPSR 

EPSI(NMH)=EPSI(NMH)-DEPSI 
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2  DO  550IX=1,NHE 
EF(NHE.EQ.  1)  GO  TO  3 
EPSR(NMH)=EPSR(NMH)+DEPSR 
EPSI(NMH)=EPSI(NMH)+DEPSI 
H(^IMH)=H(NMH)+DH 

3  WRTTE(6,621)  DC 
DO  10I=1,NMAT 

EPSC(IH>CMPLX(EPSR(I)JEPSia)) 
10  WRTrE(6,610)I,EPSR(I)£PSI(I),H(I) 
IF  f4MAT.EQ.NLAY)  WRITE(6,607) 
IF  f4MAT.GT.NLAY)  WRITE(6,608) 
FREQ=FINT-DF 
DO  500N=1,NF 
FREQ=FREQ+DF 
WRITE(6,611)FREQ 

NONDIMENSIONALIZATION 


AKP=0.53 195D0*FREQ 
BNON=B*AKP 


ANON=A*AKP 
IFfTYP.EQ.2)  GO  TO  40 
PIA=PI/ANON 

DIVR=DSQRT(1.0D0-PIA*PIA) 

10  DO  50I=1,NLAY 
50HNON(I)=H(I)*AKP 
DO  55  1=1,4 

55  ZSUM(I)=DCMPLX(0.0D0,0.0D0) 


INTEGRATION  LOOP  FORZETA 


DO  400I=1,NZET1 
AI=I 

RI=AI*ANI 

RIM=1.0D0-RI 

RKI=1.0D0/RIM 

WTI=RKI/RIM 

RKI=RKI*RI*ETZET 

RKIS=RKI*RKI 

CALCULATE  LAMBDA  ETC  FOR  EACHETA 

DO  110L=1,NMAT 

EXR=EPSR(L)-RKIS 

EXI=EPSI(L) 

C  ALLSQROOT  (EXR.EXI) 

EXRA(L)=EXR 

EXIA(L)=EXI 

AMB(L)=DCMPLX(EXR,EXI) 

AD(1  ,L)^\MB(L)/EPSC(L) 
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110  AD(2,L)=1.0AMB(L) 

LX=NLAY+1 
DO  115K=1,KMAX 

IF(NLAY.EQ.NMAT)  AD(C,LX)=OCMPLX(O.ODO,O.ODO) 
EF(NLAY.LT.NMAT)  AD^LX)=ADtfC,LX)*SMULT 
115  CONTINUE 
DO  120L=1,NLAY 
AMHI=EXIA(L),HNON(L)*2.0D0 
IF<AMHI.GT.-30.0D0)  GO  TO  118 
TN(L)=SMULT 
GO  TO  120 

1 18  AMHR=EXRA(L)*HNON(L) 

IFflTN.EQ.O)  GO  TO  1 19 
WRITE(6,650)L,AMHR,AMHI 

119  CASH=DCOS(AMHR) 

SASH=DSIN(AMHR) 

EAMI=DEXP<AMHI) 

AMP=  1 ,0-E  AMI 
AMQ=  1 ,0+E  AMI 

TOPl=CASH*AMP+SMULT*SASH*AMQ 
BOTl=SMULT*CASH*AMQ-SASH*AMP 
TN(L)=TOP  1/BOT 1 
IFflTN.EQ.O)  GO  TO  120 
WRITE(6,65 1)I,L,TN(L) 

120  CONTINUE 

DO  200K=1,KMAX 
IF<NLAY.GT.l)  GOTO  130 
TOP=AD(K,l) 

BOT=- AD(K,  1)*  AD(K,  1  )*TN(  1) 

BF(NLAY.EQ.NMAT)  GOTO  125 
TOP=TOP+AD(K,LX)*TN(l) 

BOT=BOT+AD(K,  1)*  ADfC,LX) 

125  GO  TO  190 

130  IF(NLAY.GT.2)  GO  TO  140 
TOP= AD(K,2)*(AD(K.  1 )- AD(K,2)*TN(  1  )*TN(2)) 

BOT=-AD(K,  1)*  AD(K,2)*(  AD(K,  1  )*TN(  1 )+ AD(K,2)*TN(2» 
IF(NLAY.EQ.NMAT)  GO  TO  135 

TOP=TOP+AD(K,LX)*(AD(K.  1  )*TN(2)+ AD(K,2)*TN(  1 )) 
BOT=BOT + AD(K,  1  )*  AD$C,LX)*(  AD(K.2)-AD(K.  1  )*TN(  1  )*TN(2)) 
135  GO  TO  190 

140  IF(NLAY.GT.3)  GO  TO  150 
TOP  1 = AD(K.  1  )*(AD(K.2)-AD(K.3)*TN(2)*TN(3)) 
TOP2=AD(K,2)*TN(l)*(AD(K,2)*TN(2)+AD(K,3)*TN(3)) 
TOP=(TOP  1 -TOP2)*  AD(K,3) 

BOT 1 = AD(K,  1  )*TN(  1  )*(AD(K,3)*TN(2)*TN(3  )-AD(K,2)) 
BOT2=AD(K,2)*(AD(K,2)*TN(2)+AD(K,3)*TN(3)) 
BOT=(BOTl-BOT2)*  AD(K,  1)*AD(K,3) 

IFP.AY.EQ.NMAT)  GOTO  145 
T0P1=AD(K,1)*(AD(K,2)*TN(3)+AD(K,3)*TN(2)) 
TOP2=AD(K.2)*TN(l)*(AD(K,2)*TN(2)*TN(3)-AD(K,3)) 
TOP=TOP+AD(K,LX)*(TOP  1 -TOP2) 

BOT  1  =-AD(K,  1)*TN(  1  )*(  AD(K,2)*TN(3 )+ AD(K,3  )*TN(2)) 
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BOT2=AD(K,2)*(AD(K,2)*TN(2)*TN(3)-AD(K,3)) 

BOT=BOT+(BOT  l-BOT2)*AD(K,  1)*AD(K,LX) 

145  GO  TO  190 
150  WRITE(6,612) 

STOP 

190  DDT(K)=TOP/BOT 
IFflDDT.EQ.O)  GO  TO  200 
WRITE(6,622)I,K,DDT(K) 

200  CONTINUE 
SUM1=DCMPLX(0.0D0,0.0D0) 

SUM2=SUM1 

INTEGRATION  OVER  PHI  FOHTYP=l  AND  BESSEL  FUNCTIONS  FOR  CO-AX 
FOR  EACHZETA 

DO  300J=T,MPHI1 
PHI=AMI*(J-1) 

PHSS=DSIN(PHI) 

IFJTYP.EQ.l)  GO  TO  210 
ARGl=RKI*BNON*PHSS 
ARG2=RKI*ANON*PHSS 
ARG=DCOS(ARGl)-DCOS(ARG2) 

ARGC=DCMPLX(ARG,0.0D0) 

GO  TO  250 
210  PHCC=DCOS(PHT) 

ARGl=RKI*ANON*PHCC 

ARG2=DCOS(0.5D0*ARGl)/(PI2-ARGl*ARGl) 

ARG2S=ARG2*ARG2 
IF(J.NE.l)  GO  TO  212 
ARG 1  S=0.25D0*BNON*BNON*RKI 
GO  TO  215 

212  ARG  1  =RKI*BNON*PHSS 
ARG2=DSIN(0.5D0*ARG1)/(PHSS*RKI) 

ARG  1  S=ARG2*ARG2*RKI 
215  ARG=ARG1S*ARG2S 

ARGC=(DDT(1)*PHSS*PHSS+DDT(2)*PHCC*PHCC)*ARG 
250  IF((J.EQ.  l)OR.(J.EQ.MPHIl))  ARGC=ARGC/2.0D0 
SUM  1  =SUM  1  + ARGC 
JJ=J/2 
JJ2=2*JJ 

IF(T.EQ.JJ2)  GO  TO  255 
SUM2=SUM2+2 .0D0  *  ARGC 
255  CONTINUE 
300  CONTINUE 
EFflTYP.EQ.l)  GO  TO  305 
C 

C  MODIFICATION  REQUIRED  FOR  CO-AX 
C 

XC1=SUM1*SUM1*DDT(1)/RKI 
XC2=SUM2*SUM2*DDT(  1  )/RKI 
SUM1=XC1 
SUM2=XC2 
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C  SUM  OVERZETA 
C 

305  SUM1=SUM1*WTI 
SUM2=SUM2*WTI 
ZSUM(1)=ZSUM(1)+SUM1 
ZSUM(2)=ZSUM(2)+SUM2 
E=I/2 
112=2*11 

IF<[.NE.II2)  GO  TO  400 
ZSUM(3)=ZSUM(3)+2.0D0*SUM1 
ZSUM(4)=ZSUM(4)+2.0D0*SUM2 
400  CONTINUE 
WRITE(6,613) 

DO  405 1=1,4 

ZSUM(I)=ZSUM(I)*ETZET*  ANI 
IFflTYP.EQ.l)  GO  TO  402 
ZSUM(I)=ZSUM(I)/(AM*AM) 

GO  TO  405 

402  ZSUM(I)=ZSUM(I)*AMI 
405  WRITE(6,6 1 4)I,ZSUM(I) 

C 

C  CALCULATE  ADMITTANCE,  REFLECTION  COEFFICIENT 

C 

IF<JTYP.EQ.2)  GO  TO  410 
SMY=SMULT*32.0D0*A*ZSUM(l)/p*DIVR) 
CPY=-SMY*AKP*DIVR/fREQ*3 1.917D0) 

GO  TO  420 

410  SMY=SMULT*ZSUM(  1  )/(CPSC*DLOG(B/ A)) 

CPY=SMY*  AKP*CPSC/fREQ*3 1 .9 17D0) 

420  WRITE(6,6 1 5)SMY 
WRITE(6,616CPY 
ARGC=(  1 .0D0-SMY)/(1 .0D0+SMY) 

XR1=DREAL(ARGC) 

XR2=DIMAG(ARGC) 

G  AM=DSQRT(XR  I  *XR  1  +XR2*XR2) 

C 

C  VAX  -  VAX  FORTRAN  RETURNS  DEGREESPHZ=DTAN2D(XR2,XR1) 
CPC  -MS  FORTRAN  RETURNS  RADIANSPHZ=DTAN2(XR2,XR1) 

C 

C  PHZ=DATAN2D(XR2,XR1) 

PHZ=DATAN2(XR2,XR1) 

C 

C  PC  -  CONVERT  RADIANS  TO  DEGREES  PC  VERSION  ONLY 
C 

PHZ  =  180.0*PHZ/PI 
C 

WRITE(6.617JCR1,XR2,GAM,PHZ 
IF<NF.GT.l)  WRITE(7,618FREQ,SMY,CPY 
IF(NF.EQ.l)  WRITE(7,618)HNMH),SMY,CPY 
RLS=-20.0D0*DLOG(GAM) 

WRITE(6,625]RLS 
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VSWR=20.D0*DLOG  10((1  .ODO+GAM)/(1  .ODO-GAM)) 

WRJTE(6,626)/SWR 

IF/NF.GT.  1)WRITE(8,618J:REQ,XR1,XR2,VSWR,PHZ 
IF<NF.EQ.  1  )WRITE(8,6 18)H(NMH)^CR  1  ,XR2,  VS  WR,PHZ 
500  CONTINUE 
550  CONTINUE 
CLOSE(JNIT=5) 

CLOSE(JNIT=6) 

STOP 

C 

601  FORMAT(/AX,'INPUT  DATA  FOR  CALCULATION  OF  NORMALIZED 
1' ADMITTANCE  AND  REFLECTIONCOEFFICIENTS\///,4X, 

l'NO  OF  FREQUENCIES- ,14,’  STARTINGVALUE=',E12.4,'INC=', 

1E12.4,/) 

602  FORMAT(4X, 'GUIDE  OR  CO-AX  CABLE  DIMENSIONSAB=',2E12.4,/) 

603  FORMAT/lX,'INTEGRN  PARAMETERSNZET A=',I5, 
l'MPHI=',I5,'ETZET=',E12.4,/) 

604  FORMAT^X/CASE  OF  RECTANGULAR  GUIDHTYP=',I4,/) 

605  FORMAT/IX/CASE  OF  CO-AXIAL  CABLEITYP=',I4,/) 

606  FORMAT^X,'NO  OFMTLS=',I4,’NO  OFLAYERS=',I4,/) 

607  FORMAT(///4X/LAYERS  ARE  BACKED  BY  CONDUCTING  SHEET',///) 

608  FORMAT(///)X, 'LAYERS  ARE  BACKED  BY  HALF  SPACE',///) 

609  FORMAT(4X,'DIELECTRICCONST  FOR  CABLECORE=',2E12.4,/) 

610  FORMAT^X/DIELECTRICCONST  FORMTL  ',13,'  =',2E12.4, 
lTfflCKNESS=',E12.4) 

611  FORMAT////)  OX/RESULTS  FORFRQUENCYINGHZ=',E12.4,///) 

612  FORMAT(4X,'NLAY.GT.3  EXPRESSION  NOT  AVAILABLE  IN  CODE’,/) 

613  FORMAT///4X, 'RESULTS  FOR  ACCURACY  CHECK  FOR  INTEGRALS',/, 
14X,'1-N*M,  2-N*(M/2),  3-(N/2)*M,  4-{N/2)*(M/2)',//) 

614  FORMAT/IX, 'VALUE  OF  INTEGRALNO',I3,2E12.4) 

615  FORMAT(//4X,'CALCULATIONS  BELOW  ARE  FOR  ACCURATE  VALUE', 
1  //,4X,'NORMALIZED  ADMITT  ANCE-.2E  1 2.4,/) 

616  FORMAT/IX, 'ACTUAL  ADMITTANCE=',2E12.4,/) 

617  FORM  AT/4X,'REFL  COEF=',2E  1 2.4,'MAGNITUDE=',E  12.4,/, 

14X,'PHASE  IN  DEG=',E  12.4,/) 

618  FORMAT(IX,5E12.4) 

620  FORMAT/IX,'NO  OF  PROPERTY  OR  LAYUPSETS=',I4,/) 

621  FORMAT(//,4X,'***  RESULTS  FORLAYP  SET',13,'  GIVEN  BELOW') 

622  FORMAT/lX,'IZET=',I4,'  DDT(',I2,')=',2E12.4) 

625  FORMAT(4X,'RL- .E12.4,/) 

626  FORMAT(IX,'VSWR=',E12.4,/) 

650  FORMAT  (1  X.'L  A  YER  'I2.'LMDARH.2LMDAIH=\2E  1 2.4) 

65 1  FORMATC  X,'IZET'.I4,'LAYER',I2,'FNCTN=',2E  1 2.4) 

END  OF  MAIN  PROGRAM 
END 

SUBROUTINE  SQROOT 

SUBROUTINE  SQROOT (PRE  AL,PIMAG) 

IMPLICITREAL*8  (AH,0-Z) 
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*)  w 


SRSQ=DSQRT(PREAL*PREAL+PIMAG*PIMAG) 

SR=DSQRT(SRSQ) 

THETA=-DASIN(RIMAG/SRSQ) 

IF(PREAL.LT.0.0D0jrHETA=-THETA-3. 1415926D0 

THETA=THETA/2.0D0 

PREAL=SR*DCOS(THETA) 

PIMAG=SR*DSIN(THETA) 

RETURN 

END 
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COMPARISON  WITH  BAKHTIARI'S  RESULTS  FOR  RECTANGULAR  WAVE  GUIDE 


Figures  C1-C6  show  the  variations  of  normalized  conductances  and  susceptances  for  a 
single  layer  with  the  layer  thickness  d  (backed  by  a  conductor)  for  three  values  of  the  di¬ 
electric  properties  at  a  frequency  of  10  GHz.  The  wave  guide  dimensions  are  a  =  0.9" 
and  b  =  0.4".  It  may  be  noted  that  the  present  results  are  the  same  for  series  as  well  as 
integral  solutions.  They  are  compared  with  those  reported  by  Bakhtiari,  S.  (Figures  3.3a 
and  b  in  "Open-Ended  Sensors  for  Microwave  Nondestructive  Evaluation  of  Layered  Com¬ 
posite  Media",  Ph.D.  Dissertation,  Colorado  State  University,  Fort  Collins,  CO,  1992).  It 
may  be  noted  that  although  the  trend  of  the  results  are  the  same  as  those  from  the  pres¬ 
ent  calculations,  some  of  the  peak  values  (occurring  for  some  specific  values  of  d)  reported 
by  Bakhtiari  are  lower  than  those  obtained  here.  The  likely  reason  is  that  the  integrals  (and 
the  series)  are  evaluated  in  this  study  with  a  greater  accuracy. 

Some  cases  considered  by  Bakhtiari,  et  al  (Bakhtiari,  S.,  Ganchev,  S.I.,  and  Zoughi,  R., 
"Microwave  Swept  Frequency  Optimization  for  Accurate  Thickness  on  Dielectric  Property 
Monitoring  of  Conductor-Backed  Composites",  Materials  Evaluation,  p.  740,  June  1993) 
were  also  studied  in  the  present  program.  The  variation  of  the  phase  of  the  complex  reflec¬ 
tion  coefficient  with  frequency  at  x-band  (8.2  to  1 2.4  GHz,  the  wave  guide  dimensions  are 
a =0.9",  b=0.4")  over  chosen  band  widths  were  studied  in  the  paper.  Figure  C7  (d=4.28 
mm,  £r  =  12.4-j2.4)  shows  the  results  over  9.5  to  1 1 .0  GHz,  which  are  the  same  as  those 
in  Figure  3  (measured  and  calculated)  in  the  paper.  Results  in  Figure  C8  (d  =  1.95,  2  and 
2.05  mm,  s,  =  20.0-j0.02)  are  also  close  to  those  given  in  Figure  4  of  Bakhtiari's  work. 
Figure  C9  (d  =  2.08,  2.18  mm,  £r=  1 2.4-j2.4)  compare  well  with  the  results  in  Bakhtiari's 
Figure  5.  However,  some  of  the  results  given  in  Figure  7  (d  =  2.0mm)  of  the  paper  do  not 
compare  with  the  present  calculations  shown  in  Figure  CIO.  Although  the  measured  and 
calculated  values  (in  the  paper)  for  £,  =  20-10.02  compare  extremely  well  with  those  in  Fig¬ 
ure  CIO  over  the  entire  frequency  range,  there  are  some  differences  in  the  two  results  for 
£r  =  19-j0.02.  For  sr  =  21-j0.02,  however,  the  results  (only  calculated  values  are  given  in 
the  paper)  are  drastically  different.  <j>  showing  a  drastic  change  (from  +180°  to  -180° 
with  increasing  frequency)  in  Bakhtiari's  calculations  at  a  frequency  of  about  8.4  GHz. 
This  change  is  absent  in  Figure  CIO,  which  show  smooth  variations  of  the  phase  with  fre¬ 
quency  as  well  as  the  dielectric  constant.  Calculations  were  performed  over  wider  fre- 
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I* 


quency  ranges  and  various  other  values  of  dielectric  constants  and  the  results  are  given  in 
Figures  C11-C15.  No  drastic  changes  in  <|>  were  observed  except  at  a  frequency  range  of 
16.4-17.2  GHz,  but  the  changes  are  from  -180  to  180  with  increasing  frequency  (similar 
to  that  shown  in  Figure  C7)  and  they  were  found  to  occur  for  all  three  values  of 
sI(19-j0.02,  20-j0.02  and  21-j0.02).  It  appears  that  there  may  be  some  error  in  this  result 
reported  by  Bakhtiari. 


e*  =  10.0,  fan  6  =  0.1 


Figure  Cl.  Real  Part  Normalized  Conductance  as  a  Function  of  Slab  Thickness  for  Lossy 

Dielectric  Materials 
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tr  -  10.0,  fan  6  =  0.1 


Figure  C2.  Imaginary  Part  Normalized  Susceptance  as  a  Function  of  Slab  Thickness  for 

Lossy  Dielectric  Materials 


er‘  =  5.0,  tan  <5  =  0.1 


Figure  C3.  Real  Part  Normalized  Conductance  as  a  Function  of  Slab  Thickness  for  Lossy 

Dielectric  Materials 
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tT  =  20  -  i  0.02 


Figure  C8.  Theoretical  Phase  for  Low-Loss  TBC  Sheets  with  Various  Thicknesses 


c,  =  12.4  -  i  2.4 


Figure  C9.  Theoretical  Phase  for  Low-Loss  TBC  Sheets 
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Frequency  (GHz) 


Figure  CIO.  Theoretical  Phase  for  Low-Loss  TBC  Sheets  with  Various  Permittivities 


d  =  2.0  mm 


Figure  Cl  1.  Theoretical  Phase  for  Low-Loss  TBC  Sheets  with  Various  Permittivities 


2.0  mm 


Figure  Cl  2.  Theoretical  Phase  for  Low-Loss  TBC  Sheets  with  Various  Permittivities 


d  =  2.0  mm 


Figure  Cl  3.  Theoretical  Phase  for  Low-Loss  TBC  Sheets  with  Various  Permittivities 
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d  ~  2.0  mm 


Figure  Cl  4.  Theoretical  Phase  for  Low-Loss  TBC  Sheets  with  Various  Permittivities 


d  =  2.0  mm 


Figure  Cl  5.  Theoretical  Phase  for  Low-Loss  TBC  Sheets  with  Various  Permittivities 
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RFSULTS  FOR  VARYING  AIR  GAPS  BETWEEN  A  RECTANGULAR  WAVE 
Gy|DE  AND  DIELECTRIC  SLABS  OF  VARYING  THICKNESSES  BACKED  BY  AIR 

Unavoidable  air  gaps  exist  between  the  microwave  sensors  and  the  dielectric  slab  being 
interrogated.  This  appendix  gives  the  results  of  a  systematic  study  of  the  effects  of  air 
gaps  between  a  rectangular  wave  guide  and  the  slab  thickness  d  for  ^=2.8  and  tan5  =  .05. 
The  slab  is  backed  by  air.  The  results  (actual  admittances  in  ohm'1)  are  given  in  the  fol- 


lowing  figures. 

Figure 

Plots 

d,  inch 

D1 

Conductance  vs.  air  gap  and  frequency 

.  0.075" 

D2 

Susceptance  vs.  air  gap  and  frequency 

0.075" 

D3 

Conductance  vs.  air  gap  and  frequency 

0.125" 

D4 

Susceptance  vs.  air  gap  and  frequency 

0.125" 

D5 

Conductance  vs.  air  gap  and  frequency 

0.235" 

D6 

Susceptance  vs.  air  gap  and  frequency 

0.235" 

D7 

Conductance  vs.  air  gap  and  frequency 

0.38" 

D8 

Susceptance  vs.  air  gap  and  frequency 

0.38" 

D9 

Conductance  vs.  air  gap  and  frequency 

Half  Space 

D10 

Susceptance  vs.  air  gap  and  frequency 

Half  Space 

Dll 

Susceptance  vs.  conductance  for  varying  air  gaps  and  frequencies 

0.075" 

D12 

Susceptance  vs.  conductance  for  varying  air  gaps  and  frequencies 

0.125" 

D13 

Susceptance  vs.  conductance  for  varying  air  gaps  and  frequencies 

0.235" 

D14 

Susceptance  vs.  conductance  for  varying  air  gaps  and  frequencies 

0.38" 

D15 

Susceptance  vs.  conductance  for  varying  air  gaps  and  frequencies 

Half  Space 

The 

results  show  the  strong  influence  of  the  air  gap  thickness  on  the 

susceptance 

when  the  thickness  is  small  (<0.05").  On  the  other  hand,  conductance  of  slabs  of  small 
thickness  (<0.125")  is  not  that  sensitive  to  small  air  gap  thicknesses.  For  slabs  of  larger 
thicknesses  (>0.235"),  the  conductance  is  also  sensitive  to  small  air  gaps.  Frequency  in¬ 
fluences  the  results  most  for  certain  ranges  of  air  gaps  (usually  for  small  or  large  values). 
In  general,  it  appears  that  smallest  possible  gap  is  desirable  for  obtaining  better  resolutions 
for  imaging  purposes,  but  slightly  higher  values  (<0.02")  may  be  acceptable. 
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Figure  D1 .  Conductance  of  a  Dielectric  Slab  of  Thickness  0.075  inch  and 
er  =2.8,  tan8  =  .05  Backed  by  Air  for  Varying  Air  Gap 


Figure  D2.  Susceptance  of  a  Dielectric  Slab  of  Thickness  0.075  inch  and 
Er  =2.8,  tan5  =  .05  Backed  by  Air  for  Varying  Air  Gap 
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Figure  D5.  Conductance  of  a  Dielectric  Slab  of  Thickness  0.235  inch  and 
s,  =2.8,  tanS  =  .05  Backed  by  Air  for  Varying  Air  Gap 
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Figure  D6^  Susceptance  of  a  Dielectric  Slab  of  Thickness  0.235  inch  and 
s,  =2.8,  tan5  =  .05  Backed  by  Air  for  Varying  Air  Gap 
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Figure  D7.  Conductance  of  a  Dielectric  Slab  of  Thickness  0.380  inch  and 
Sr  =2.8,  tan5  =  .05  Backed  by  Air  for  Varying  Air  Gap 


Figure  D8.  Susceptance  of  a  Dielectric  Slab  of  Thickness  0.380  inch  and 
er  =2.8,  tan5  =  .05  Backed  by  Air  for  Varying  Air  Gap 
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Figure  Dll.  Relationship  of  Susceptance,  Conductance,  &  Air  Gap  (d)  on  a  Dielectric  Slab 
of  Thickness  0.075  inch  and  s,  =2.8,  tan8  =  .05  Backed  by  Air 
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Figure  D12.  Relationship  of  Susceptance,  Conductance,  &  Air  Gap  (d)  on  a  Dielectric  Slab 
of  Thickness  0.125  inch  and  sr  =2.8,  tan5  =  .05  Backed  by  Air 
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Figure  D13.  Relationship  of  Susceptance,  Conductance,  &  Air  Gap  (d)  on  a  Dielectric  Slab 
of  Thickness  0.235  inch  and  Sr  =2.8,  tan8  =  .05  Backed  by  Air 


Figure  D14.  Relationship  of  Susceptance,  Conductance,  &  Air  Gap  (d)  on  a  Dielectric  Slab 
of  Thickness  0.380  inch  and  sr  =2.8,  tand  =  .05  Backed  by  Air 
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Figure  D15.  Relationship  of  Susceptance,  Conductance,  &  Air  Gap  (d)  on  a  Dielectric 

Half  Space,  sr  =2.8,  tan8  =  .05 
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INTRODUCTION 


Dielectric  properties  of  layered  materials  are  often  influenced  by  air  gaps  or  delamina¬ 
tions,  which  may  be  present  between  the  layers.  Use  of  microwave  NDE  for  detecting  the 
air  gap  at  such  delaminations  has  been  suggested  by  some  investigators.  In  this  approach 
the  sensor  (a  rectangular  wave  guide)  is  placed  in  front  of  a  coupling  dielectric  sheet, 
which  is  backed  by  the  specimen  being  interrogated.  An  unintentional  air  gap  is  usually 
present  between  the  coupling  material  and  the  specimen  surface.  This  is  called  the  lift-off 
effect  and  one  first  attempts  to  find  the  direction  of  the  change  in  the  tip  of  the  vector 
(complex  reflection  coefficient)  due  to  small  lift-off  (air  gap).  This  direction  (that  of  the 
change)  is  called  the  lift-off  direction.  For  the  purpose  of  detecting  the  presence  or  ab¬ 
sence  of  a  delamination  type  air  gap  inside  the  specimen,  the  change  in  the  component  of 
the  complex  reflection  coefficient  vector,  which  is  orthogonal  to  the  lift-off  direction,  is 
normally  used,  so  as  to  minimize  the  lift-off  effects  produced  by  wobbling  of  the  wave 
guide  and  the  coupling  sheet  over  the  specimen  during  scanning.  At  certain  frequencies 
the  lift-off  is  a  minimum  and  the  orthogonal  component  is  a  maximum.  This  Appendix  de¬ 
scribes  the  results  of  a  parametric  study  performed  to  find  the  frequency  ranges  which 
should  be  used  so  as  to  minimize  the  lift-off  effect  in  support  of  microwave  NDE  studies  at 
NSWC. 
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LIFT-OFF  DIRECTION 


The  lift-off  direction  is  determined  by  considering  the  following  cases,  i.e.,  the  presence 
of  a  small  air  gap  c  of  varying  thicknesses  and  no  air  gap  between  the  coupling  material  B 
(next  to  the  rectangular  wave  guide)  and  the  sample  material  D  as  shown  in  Figure  El.  For 
this  study,  the  following  parameters  were  chosen. 

Region  A  -  x-band  wave  guide  a  =  0.9",  b=0.4" 

Region  B  -  Coupling  dielectric  sheet,  thickness  =  0.375",  e,=2.7,  tan5=0.05 
Region  C  -  Unintentional  air  gap  of  thickness  t  varied  from  0  to  0.02" 

Region  D  -  Sample,  £,  =  4.5,  tanS=0.05,  thickness  =  4" 


C 


_l 
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e> 

3> 

Ai<*~ 

1 

Figure  El.  Geometry  to  Determine  the  Lift-Off  Direction 
A  -  Wave  Guide 
B  -  Coupling  Material 
C  -  Air  Gap 
D  -  Sample 

The  objective  was  to  find  the  trajectory  of  the  tip  of  the  reflection  coefficient  factor  (X,Y,  R 
=  X+iY)  as  the  unintentional  air  gap  thickness  was  changed  (using  the  code  ADMIT  for  cal¬ 
culations)  and  the  tangent  to  this  trajectory  for  zero  air  gap  thickness  (t  =  0).  It  was  found 
that  a  line  drawn  from  the  tip  for  zero  air  gap  thickness  to  that  for  a  gap  thickness  of 
t  =  0.01 "  is  adequate  to  define  the  tangent  or  the  lift-off  direction.  Let 
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RXX,  RYY  =  X.Y  for  t  =  0 
RXZ,  RYZ  =  X,Y  for  t  =  0.0 1 " 

DX  =  RXZ-RXX 
DY  =  RYZ-RYY 

The  direction  cosines  of  the  lift-off  direction  is 

DXX  =  DX/DL  <E'2> 

DYY  =  DY/DL 
DL  =  (DX2  +  DY2),Q 

The  direction  cosines  of  the  orthogonal  vector  (to  lift-off  vector)  is 

PDX  =  -DYY  (E'3> 

PDY  =  DXX 

The  projection  of  the  reflection  coefficient  vector  for  t=0  on  to  the  lift-off  direction  is 

RDXZ  =  RXX-DXX  +  RYY-DYY  (E-4) 

and  on  to  the  orthogonal  direction  is 

RDPZ  =  RXX-PDX  +  RYY*PDY  (E-5) 

Note  that  all  the  quantities  defined  above  are  functions  of  frequency. 
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MINIMIZATION  OF  LIFT-OFF  EFFECT 


This  minimization  is  desirable  to  increase  the  accuracy  of  prediction  of  presence  of  de¬ 
laminations  within  the  sample  D  at  different  locations  (Figure  E2).  The  following  cases 
were  considered  as  shown  in  Figure  E2,  but  the  thickness  of  the  unintentional  air  gap  C 
was  chosen  as  zero. 

Distance  of  air  gap  in  D  from  C,  x  =  1 ",  2",  3" 

Air  gap  thickness  8  =  0.005",  0.01 ",  0.02" 


Figure  E2.  Geometry  for  Study 
A  -  Wave  Guide 
B  -  Coupling  Material 
C  -  Unintentional  Air  Gap  is  Absent 
D  -  Sample  with  Delamination 

The  reflection  coefficient  vectors  (X,  Y)  for  each  of  these  cases  were  computed  using  the 
code  ADMIT  and  their  projections  on  lift-off  and  orthogonal  directions  were  calculated. 

Lift-Off:  DRDX  =X*DXX  +  Y»DYY  (6) 

Orthogonal:  DRPD=X»PDX  +  Y»PDY 

Finally,  changes  of  these  projections  from  those  for  no  defect  case  (equation  (E4)  and  (E5)) 
were  determined. 

Lift-Off:  DRLF=DRDX  -  RDXZ 
Orthogonal:  DRPF=DRPD  -  RDPZ 
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RESULTS 


Values  of  DRLF  and  DRPF  for  5  =  0.005",  0.01",  and  0.02"  are  plotted  against 
frequency  (9  to  1 1  GHz)  in  Figures  E3,  E4,  and  E5  for  the  defect  locations  x  =  1",  2", 
and  3",  respectively.  It  may  be  noted  that  at  the  frequencies  the  orthogonal  compo¬ 
nent  DRPF  are  maximum,  the  lift-off  components  DRLF  are  minimum.  As  the  distance 
x  of  the  defect  from  the  coupling  material  is  increased,  the  components  reduce  in  mag¬ 
nitude  (an  expected  effect)  and  the  maximum  values  are  attained  at  an  increasing 
number  of  frequencies.  The  frequencies  at  which  the  orthogonal  components  are 
maximum  should  be  chosen  for  interrogation  to  determine  existence  of  possible  delami- 
hations  in  the  sample.  Figure  E6  shows  the  variations  of  the  projections  of  the  reflec¬ 
tion  coefficient  vector  (RDXZ,  RDPZ)  for  no  defect  case  on  to  the  lift-off  and  orthogonal 
directions  for  various  frequencies. 


Figure  E3.  Frequency  vs.  Lift-Off  and  Orthogonal  Components 
Air  Gap  at  1 "  from  Front  Side  in  Material  2 
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Lift-Off  &  Orthogonal  Component 


Figure  E4.  Frequency  vs.  Lift-Off  and  Orthogonal  Components 
Air  Gap  at  2"  from  Front  Side  in  Material  2 


Figure  E5.  Frequency  vs.  Lift-Off  and  Orthogonal  Components 
Air  Gap  at  3"  from  Front  Side  in  Material  2 
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Figure  E6.  Projection  of  Reflection  Coefficient  for  Zero  Defect  on  to  Lift-Off 

and  Orthogonal  Direction 
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INTRODUCTION 

This  appendix  gives  some  comparison  of  results  obtained  using  the  computer  code 
ADMIT  for  the  co-axial  line  with  the  results  presented  by  Bakhtiari,  et  al.,  and  Li,  et  al. 
The  results  show  excellent  comparison  and  the  fact  that  the  computer  code  is  free  from 
any  bugs  and  the  results  can  be  used  with  confidence  for  selecting  frequencies  and  other 
parameters  for  interrogating  specimens. 

Figures  FI  -  F10  compare  the  results  with  those  reported  by  Bakhtiari,  et  al.  (Bakhtiari, 
S.,  Ganchev,  S.  L,  and  Zoughi,  R.,  "Analysis  of  Radiation  from  an  Open-Ended  Co-Axial 
Line  into  Stratified  Dielectrics",  IEEE  Transactions  on  Microwave  Theory  and  Techniques, 
Vol.  42,  No.  7.  July  1994).  There  are  no  differences  in  the  results. 

Figure  F1 1  compares  the  results  with  those  reported  by  Li,  et  al.,  (Li,  C.L.,  and  Chen, 
K.M.,  "Determination  of  Electromagnetic  Properties  of  Materials  Using  Flanged  Open-Ended 
Co-Axial  Probe  -  Full  Wave  Analysis",  IEEE  Transactions  on  Instrumentation  and  Measure¬ 
ment,  Vol.  44,  No.  1,  p.  19,  February  1995).  The  results  are  found  to  differ  at  higher  fre¬ 
quencies,  either  due  to  the  effect  of  small  damping  included  in  the  present  analysis  or  the 
consideration  of  higher  modes  considered  by  Li,  et  al.  Most  likely,  the  later  effect  plays  a 
more  dominant  role. 
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Figure  FI 


u 


f 


=  10,  c'rt  =  1  (air),  f  =  5  GHz,  a  =  1.18  mm,  b  =  3.62  mm 


Return  Loss  Versus  Normalized  Slab  Thickness  for  a  Dielectric  Layer 
Backed  by  Free-Space 


e’,,  =  10,  e*rt  =  1  (air),  f  =  5  GHz,  a  =  1.18  mm.  b  =  3.62  mm 


Figure  F2.  Phase  Versus  Normalized  Slab  Thickness  for  a  Dielectric  Layer 

Backed  by  Free-Space 


(deg) 


Figure  F3.  Return  Loss  Versus  Normalized  Slab  Thickness  for  a  Dielectric  Layer 

Backed  by  a  Conductor 


Figure  F4.  Phase  Versus  Normalized  Slab  Thickness  for  a  Dielectric  Layer 

Backed  by  a  Conductor 
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=  10.  e’„  =  2.07,  f  =  5  GHz,  o  =  1.18  mm,  b  =  3.62  mm 


0.25 


.y 


tr  4  +*■*’*' 


2*5^* 


Frequency  f  (GHz) 


Figure  F5.  Return  Loss  Versus  Frequency  for  a  Dielectric  Layer  Backed  by  Free-Space 


=  10,  £’„  =  2.07,  f  =  5  GHz.  a  =  1.18  mm,  b  =  3.62  mm 
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Figure  F6.  Phase  Versus  Frequency  for  a  Dielectric  Layer  Backed  by  Free-Space 
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e%,  =  10,  t'„  =  2.07,  f  =  5  GHz,  o  =  1.18  mm.  b  =  3.62  mm 


Frequency  f  (GHz) 


Figure  F7.  Return  Loss  Versus  Frequency  for  a  Dielectric  Layer  Backed  by  a  Conductor 

c’n  =  10,  c\  =  2.07,  f  ss  5  GHz,  a  =  1.18  mm,  b  =  3.62  mm 


Figure  F8.  Phase  Versus  Frequency  for  a  Dielectric  Layer  Backed  by  a  Conductor 
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c’,,  =  10.  »an<5,  =  0.01,  e'*  =  1  (air),  d,  =  0.5  mm,  a  =  1.18  mm,  b  -  3.62  mm 


Figure  F9.  Return  Loss  Versus  Frequency  for  Various  Disbonds  in  a  Two  Layer  Medium 

Backed  by  a  Conducting  Sheet 


c*rt  =  10,  tanfi,  =  0.01,  e'rt  =  1  (air),  d,  =  0.5  mm.  a  =  1.18  mm,  b  =  3.62  mm 


Figure  F10.  Phase  Versus  Frequency  for  Various  Disbonds  in  a  Two  Layer  Medium  Backed 

by  a  Conducting  Sheet 
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Reflected  Coefficient,  R  (dB) 


c’rt  =  2.3,  a  =  1.4364  mm,  b  =  4.725  mm 


Figure  F1 1 .  Reflection  Coefficient  and  Phase  of  a  Coaxial  Probe  Placed  Against  an  Infinite 

Homogeneous  Medium 
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In  this  appendix,  results  obtained  by  Liu  ("Electrical  Characterization  of  Fiber  Reinforced 
Composites  for  Mast  and  Radome  Applications,  CARDIVNSWC-TR-95/208)  were  corre¬ 
lated  with  predictions  from  the  code  ADMIT  for  investigating  the  effects  of  the  following 
factors  on  the  measurement  of  dielectric  constants  using  the  co-axial  line. 

i.  Effect  of  thickness  of  the  dielectric  layer 

ii.  Effect  of  surface  roughness 

For  low  frequencies  the  admittance  Y  of  the  probe  can  be  expressed  in  terms  of  two 
constants  Co  and  C,  which  depend  on  the  dimensions  of  the  probe.  (See  Misra,  D.,  et  al., 
IEEE  Transactions  on  Microwave  Theory  and  Techniques,  Vol.  38,  No.  1 ,  January  1 990.) 

Y  =  C0+ jaC^e' -js”)  (G1) 

«®C,(£*  +  js ') 

where  s'  is  the  dielectric  constant  of  the  material  being  interrogated  and 
e"  =  s' tan  S,  tan£  being  the  loss  tangent.  Equation  (G2)  is  often  used  in  practice,  since 
Co  is  usually  small.  The  computer  code  was  used  to  calculate  the  admittance  Y  for  the 
following  probe  with  Teflon  as  the  core  material  (sc  =  2.07)  placed  against  semi-infinite 
media  with  varying  dielectric  constants. 

a  =  inner  radius  =  0.5959  mm  =  0.02323"  (G3) 

b  =  outer  radius  =  1.81  mm  =  0.07126"  (G4) 

Results  are  given  in  Tables  G1  and  G2  for  various  frequencies  and  dielectric  constant  s  of 
the  medium.  Two  values  of  tan8  are  considered  (tan8  =  0.02  in  Table  G1  and  0.2  in  Table 
G2).  The  results  show  that  equation  (G2)  with  C|  =  2.062  can  yield  good  estimates  of 
the  dielectric  constant  s'  (with  10%  accuracy)  for  both  values  of  tan8  and  a  wide  fre¬ 
quency  range.  However,  use  of  the  same  equation  will  not  yield  a  reliable  estimate  of  s" , 
since  the  linear  relationship  with  o>  does  not  hold  for  higher  values  of  co  especially  when 
tan8  is  small.  Therefore,  the  calibration  procedure  for  evaluation  of  tan8  is  possibly  more 
complicated.  However,  since  the  objective  of  this  study  is  to  evaluate  the  effect  of  the 
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thickness  of  medium  and  surface  roughness  on  the  dielectric  constant  s' ,  we  use  equation 
(G2)  as  a  basis  for  comparison. 

EFFECT  OF  LAYER  THICKNESS 

The  loss  tangents  of  the  composites  studied  by  Liu  are  small  and  the  dielectric  con¬ 
stant  s'  are  also  not  large  (of  the  order  of  2.6  to  5.1).  The  admittances  Y  were  calculated 
for  various  layer  thicknesses  and  the  following  chosen  values  of  s'  and  tan8. 

s'  =  4.7  <G5) 

tan5  =  0.02  .  <G6> 

The  results  are  given  in  Table  G3  and  they  show  that  variations  of  layer  thicknesses  con¬ 
sidered  (Liu  used  1  to  3  layer  plates,  a  single  layer  being  0.125"  thick)  will  not  yield  any 
significant  error  in  evaluation  of  the  dielectric  constant  s' .  Estimates  of  tan8  may,  how¬ 
ever,  be  erroneous  for  high  frequencies  (©>7  GHz). 

EFFECTS  OF  SURFACE  ROUGHNESS 

To  investigate  the  effects  of  surface  roughness,  the  admittance  was  calculated  for  air 
gap  thicknesses  of  0  to  0.0024"  between  the  probe  and  the  dielectric  layer  of  thickness 
0.125".  The  properties  of  the  layer  were  chosen  to  be  the  same  as  in  equations  (G5)  and 
(G6).  Results  are  shown  in  Table  G4.  The  value  of  Yi,  reduces  by  about  12.5%  (for  ©  be¬ 
tween  1  and  11  GHz)  as  an  air  gap  of  thickness  0.0008"  is  introduced.  The  rate  of  the 
reduction  gradually  reduces  as  the  air  gap  thickness  increases  and  a  total  reduction  of 
28%  occurs  for  an  air  gap  of  0.0024".  This  change  is  in  agreement  with  that  reported  by 
Liu  for  the  material  designated  as  7781 /510A.  The  difference  between  the  measured  val¬ 
ues  of  s'  from  polished  and  unpolished  surfaces  was  found  to  be  about  35%.  The  unpol¬ 
ished  surface  (examined  by  profilometry)  show  large  undulations  and  smaller  roughness 
features  (ripples)  with  roughness  values  of  the  order  of  45  and  20  microns.  For  an  as¬ 
sumed  sinusoidal  variation,  a  RMS  value  of  20  microns  corresponds  to  an  air  gap  of 
0.0024",  the  maximum  value  of  the  gap  studied  in  Table  G4.  The  effect  on  the  real  part 
of  the  admittance  (and  s"  or  tan8)  is  more  pronounced. 
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We  also  attempted  to  Investigate  the  effect  of  an  air  gap  between  two  layers  In  addi¬ 
tion  to  an  air  gap  of  0.0024"  between  the  probe  and  the  top  layer.  Y  is  affected  margin¬ 
ally  due  to  the  second  air  gap. 
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Table  G1 .  Calculated  Values  of  Ct,  for  Various  ©  and  s' ,  tan8  =  0.02 


C,  ~  104F,  /(e'o)) 


© 

GHz 

■MCgH 

mam 

hb 

H 

Ha 

1 

2.061 

2.063 

2.063 

2.070 

2.079 

3 

2.067 

2.078 

2.084 

2.136 

2.205 

5 

2.083 

2.112 

2.126 

2.253 

2.381 

7 

2.102 

2.157 

2.182 

2.376 

2.485 

9 

2.130 

2.212 

2.248 

2.466 

2.423 

11 

2.161 

2.272 

2.318 

2.484 

2.159' 

C,  «  104Fr  l{e'<o) 


© 

GHz 

HnnH 

HB 

HB 

warn 

Ha 

H 

1 

2.064 

2.070 

2.076 

2.124 

2.225 

3 

2.114 

2.237 

2.312 

3.411 

5.633 

5 

2.285 

2.802 

3.117 

7.695 

1 6.378 

7 

2.634 

3.971 

4.784 

16.003 

35.000 

9 

3.225 

5.921 

7.540 

28.361 

58.220 

11 

4.109 

8.781 

11.525 

43.770 

80.420 
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Table  G2.  Calculated  Values  of  C,  for  Various  ©  and  s' ,  tan8  =  0.2 


C,  « 104r,  l{e'(o) 


© 

GHz 

H 

HB 

■m 

■Hi 

HI 

1 

2.061 

2.062 

2.063 

2.069 

2.077 

3 

2.067 

2.077 

2.081 

2.123 

2.172 

5 

2.081 

2.104 

2.116 

2.202 

2.260 

7 

2.097 

2.138 

2.156 

2.257 

2.236 

9 

2.118 

2.176 

2.198 

2.260 

2.052 

11 

2.142 

2.212 

2.235 

2.185 

1.719 

C,  « 1047R  /(s'co) 


© 

GHz 

s'  =2 

HB 

mpmi 

H9 

H 

m a 

1 

2.062 

2.065 

2.067 

2.083 

2.108 

3 

2.079 

2.111 

2.127 

2.324 

2.647 

5 

2.123 

2.223 

2.278 

2.916 

3.908 

7 

2.191 

2.410 

2.529 

3.861 

5.704 

9 

2.294 

2.686 

2.895 

5.100 

7.593 

11 

2.432 

3.049 

3.373 

6.456 

9.090 
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Table  G3.  Calculated  Values  of  Ci  for  Varying  Layer  Thickness  h 
e' ,  =4.7,  tan5  =  0.02 


C,  * lO4^  /(g'co) 


CD 

GHz 

ISB 

mm 

h  =oo 

1 

1.980 

2.048 

2.099 

2.063 

3 

1.988 

2.060 

2.075 

2.078 

5 

2.009 

2.090  ] 

2.114 

2.112 

7 

2.038 

2.137 

2.165 

2.157 

9 

2.082 

2.213 

2.217 

2.212  • 

11 

2.139 

2.309 

2.267 

2.272 

C,  «10 4YR/(e'<a) 

CD 

GHz 

mm 

IBH 

IH 

h=co 

1 

1.948 

2.044 

2.062 

2.070 

3 

1.979 

2.072 

2.105 

2.237 

5 

2.022 

2.159 

2.470 

2.802 

7 

2.206 

2.465 

3.951 

3.971 

9 

2.281 

2.565 

6.300 

5.921 

11 

2.760 

6.656 

9.204 

8.781 
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0.962  x  1 0'3 


0.291  x  10- 


0.491  x  10*2 


0.703  x  10- 


0.842  x  1 0'3 


0.254  x  10 


0.430  x  10‘2 


0.615  x  10 


9 

0.936  x  10-2 

0.819  x  10  2 

11 

0.119  x  lO'1 

0.104  x  101 

0.754  x  10  3 


0.228  x  10 


0.385  x  1 0'2 


0.552  x  10- 


0.734  x  10 


0.935  x  10- 


0.687  x  IQ3 


0.207  x  10 


0.350  x  10- 


0.502  x  1 0'2 


0.668  x  10 


0.849  x  10'2 


g  =  0 


0.192  x  10 


0.584  x  10-4 


0.102  x  10  3 


0.162  x  10'3 


g  =  0.0008" 


0.146  x  10 


0.444  x  10-4 


0.774  x  10  4 


0.127  x  IQ3 


g  =  0.0016" 


0.116  x  10 


0.354  x  IQ4 


0.621  x  10-4 


0.105  x  10  3 


g  =  0.0024" 


0.952  x  10 


0.291  x  IQ*4 


0.513  x  10-4 


0.889  x  10  4 


9 

0.302  x  10'3 

0.251  x  10  3 

0.217  x  10’3 

0.193  x  10  3 

11 

0.688  x  1 0  3 

0.600  x  10‘3 

0.536  x  lO'3 

0.487  x  lO3 

